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■ We present results from an extensive set of one- and two-dimensional radiation-hydrodynamic simulations 

\f-^ \ of the supernova core collapse, bounce, and postbounce phases, and focus on the protoneutron star (PNS) 

. spin periods and rotational profiles as a function of initial iron core angular velocity, degree of differential 

rotation, and progenitor mass. For the models considered, we find a roughly linear mapping between initial 
iron core rotation rate and PNS spin. The results indicate that the magnitude of the precollapse iron core 
angular velocities is the single most important factor in determining the PNS spin. Differences in progenitor 
O ' mass and degree of differential rotation lead only to small variations in the PNS rotational period and profile. 

\ Based on our calculated PNS spins, at ^200-300 milliseconds after bounce, and assuming angular momentum 

C/3 i conservation, we estimate final neutron star rotation periods. We find periods of one millisecond and shorter 

' for initial central iron core periods of below ~10 s. This is appreciably shorter than what previous studies have 

predicted and is in disagreement with current observational data from pulsar astronomy. After considering 
possible spindown mechanisms that could lead to longer periods we conclude that there is no mechanism that 
/K^ ' can robustly spin down a neutron star from ~1 ms periods to the "injection" periods of tens to hundreds of 

H I milliseconds observed for young pulsars. Our results indicate that, given current knowledge of the limitations 

of neutron star spindown mechanisms, precollapse iron cores must rotate with periods around 50-100 seconds 
to form neutron stars with periods generically near those inferred for the radio pulsar population. 

Subject headings: Stars: Rotation, Stars: Evolution, Stars: Neutron, Pulsars: General, Stars: Supernovae: 
General, Hydrodynamics, Neutrinos, Radiative Transfer 



1. INTRODUCTION 

Neutron stars are formed in the core collapse of dying mas- 
sive stars or (rarely) in the accretion-induced collapse (AIC) 
of O-Ne-Mg white dwarfs in binary systems. In both types of 
collapse events an electron-degenerate core is pushed over its 
effective Chandrasekhar mass. The collapse is halted only 
when the core density reaches ~2 x lO''^ g cm"^ and nu- 
clei undergo a phase transition to free nucleons. The equa- 
tion of state (EOS) abruptly stiffens at that point and the 
core bounces, launching a hydrodynamic shock wave into 
the supersonically infalling outer core and leaving behind a 
hot and extended protoneutron star (PNS). The bounce shock 
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Stalls almost immediately and turns into an accretion shock 
due to a combination of neutrino losses, the ram pressure 
of the infalling mantle, and the dissociation of heavy nu- 
clei into neutrons, protons and a particles. The subsequent 
evolution of the system to a potentially successful super- 
nova explosion has been under detailed investigation for many 
decades now. The fav ored "delayed neutrino mechanism" 
jBethe & Wilso nl Fml relies on shock revival due to suffi- 
cient energy deposition by neutrinos in the so-called "gain 
region" behind the stalled shock. However, using numer- 
ical simulations that incorporate Boltzmann neutrino trans- 
fer and the best available physics, a number of groups have 
shown tha t this mech anism does not succeed in spherical 
symmetrv (iMezza cappa et al..200L Liebendo r feret al, 20011 
iRampp & .Tankal l2002t iThompson et alJ l2003h . In spheri- 
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cal symmetry, the bounce shock is not revived, though an 
increase of only ~25% in the net energy deposition (all 
else being equal) would be sufficient to drive a success- 
ful explosion (jXhompsonet al. 2005). H owever, recent two- 
dimensional simulations ( Herant et al. 1994'; "Burrows et al 
p95, Frver & Hepi .2000: Janka et al.. 200.3: .Buras et al 
l200ilWa der et alJ2005h and the first three-dimensional sim- 
ulations ( iFrver & WarrMil l2004t iJanka et al.1 12005ft indicate 
that indeed multi-dimensional phenomena such as postbounce 
convection and/or rotation might lead to the conditions needed 
for a successful explosion. In the same context, Thomp- 
son et al. (2005) performed simulations in spherical symme- 
try, including a ID prescription for rotation, angular momen- 
tum transport, and viscous heating. A number of their more 
rapidly rotating models exploded due to the additional energy 
deposited behind the stalled shock by viscous dissipation. 

It is, however, possible that only stars of masses less than 
^25-3 Mf7) explode as superno vae, leaving neutron stars 
llFWe3fT999t iHegeret al.ll20()3ah . The fate of more mas- 
sive progenitor stars is more unc ertain and may be as sociated 
with under-energetic supernovae (Nomoto et al.'2 00'5l) . hyper- 
nova e, gam ma-r ay bursts (GR Bs), and/or black hole forma- 
tion (|Frver 199?; Frver & Kaloaera 2001 ; Heger et al.l2003ah . 
If shock revival and a strong supernova explosion result, no 
critical amount of mass is accreted onto the PNS. The PNS 
cools and contracts quasi-hydrostatically to its final radius of 
~10 - 15 km on a timescale of tens of seconds, radiating 
away ^ 3 x lO**^^ erg s in neutrinos feurrows & LattimeJ 19861 
iKeil & .Tanka.l995:.Pons et al.ll99^ 

Long after the supernova explosion, the still young neu- 
tron star becomes visible to obs ervers, possibly as a can onical 
rotation-powered radio pulsar ( iKaspi & Helfandl|2002ft . Ob- 
servationally, it is often difficult to directly associate pulsars 
with supernova remnants (SNRs). This is partly due to the dif- 
ference between the lifeti mes of a typical pulsar (10^~^ y ears) 
and a SNR (lO^-^ years) ( iLvne & Graham-Smithlf!998l) . In 
addition, pulsars have observed average space velocities of 
200 - 500 km s~' ( iLvne & LorimeJ 1994ft . most probably im- 
parted to them during the supernova. These "pulsar kicks" 
may be a strong indication of the multi-dimensional nature of 
the su pernova mechanism jBurrows & Havesll99^lLai et aP 
l20()!llFrveil2004tlScheck et al.l2004ft . 

Observations of radio pulsars imply that the initial rota- 
tion rate of neutron st ars is in the vicinity of tens to hun- 
dreds of milliseconds ( iKaspi & Helf and^'2002ft. The fastest 
known young pulsar is PSR J0537-6910 which is associated 
with SNR N157B and has a spin period of 16 ms and an age of 
~5 kyr If one assumes a pulsar breaking index, n, less than 
or equal to 2 an initial pulsar period < 10 ms is obtained 
( iMarshall et alj[T998l) . For the Crab pulsar (PSR B0531H-21, 
P = 33 ms), n has been measured and because the year of 
the Crab's birth is known independently (1054 A.D.), a solid 
initial spin determinatio n of Pinidai = 19 ms has been made 
jKaspi & HelfanJ2002ft . 

The empirically determined initial spin periods of pul- 
sars have to be matched to the theoretical models of stellar 
evolution, core collapse, supernova explosion, and neutron 
star cooling. Massive stars are suspected to rotate rigidly 
and rapidly during hydrogen burning on the main sequence. 

' A breaking index of 3 liolds for pure magnetic dipole radia tion and n < 2- 
3 is common for young pulsars ILvne & Graham-Smithll993) . 



Their equatorial velocities are typically ^200 km s which 
corresponds to ~10% of their Keplerian break up speed 
( lFukudalll98l . Since there is no detailed empirical knowl- 
edge of the evolution of the angular momentum distribu- 
tion inside an evolving star and because consistent multi- 
dimensional stellar evolutionary calculations have yet to be 
performed, stellar evolution theorists must rely on parameter- 
dependent, semi-phenomenological prescriptions to follow 
the angular momentum evolut ion of mass i ve sta rs in their 
spherically symmetric models. 'Heger et al.' (2000*) have pre- 
sented the first rotating presupernova stellar models that in- 
cluded a ID prescription for angular momentum transport 
and centrifugal effects (the latter only up to the end of 
core carbon burning). These authors estimated the neutron 
star spin rate by assuming that the total angular momen- 
tum contained in the precollapse iron core is conserved dur- 
ing collapse and supernova phases and is deposited com- 
pletely in a rigidly rotating neutron star of radius and mo- 
ment of inertia (I) of ~12 km a nd 0.35 Mns Rns' respec- 
tively ( iLattimer & Pra kashf'2001'). In this simple way, they 
estimated an initial neutron star period of 1 ms for their fidu- 
cial 20 Mq (ZAMS) stellar model. Given the observational 
data, this initial period is at least a factor of ten too short 
for the generic pulsar Processes that could lead to an early 
spindown of a young (proto)neutron star include: angular 
momentum redistribution by global hydrodynamic rotational 
instabilities JCentrella et alJ l2001 ; Ott et al. 2005), r-modes 
and gravit ational radiation back-reaction (iHo & Laii i2000l 
pndblom e tani200!l iWatts & Anderssoril2002t lArras et alJ 
12003), rotation-powered explosions (lAkiv amae^Ln2003F 
Akivama & Wheeler 2005), viscous angular momentum 
transport due to convection, neutrino viscosity or dis- 
sipati on of shear e n ergy stored in differential rota- 
tion (iTh ompson et alJ I2005[). magneto-c entrifugal winds 
(|Thompson''2003a'; "Thompson et al."2004ft, early magnetic- 
dipole radi ation (Ostriker & Gunn 1969), late-time fal l-back 
(iHeger et al.ll2005ft . and anisotropic neutrino emission dJankal 
l2004r 

It is also possible that slowly spinning cores are the natu- 
ral end-product of stellar evolution. Spruit & Phinnev ( 199^ 
included angular momentum transport via magnetic processes 
that yielded such large neutron star periods (P ~100 s) that the 
authors had to rely on subsequent spin-up by off-center birth 
kicks to obtain spin periods matching the observations. With 
a similar approach for magnetic angular momentum transport 
during stellar evolution and the simple neut ron star spin es- 
timate described above. iHeger et al] i2003bft and lHeger et alJ 
( 12005ft obtained neutron star periods of ~10 ms. However, 
given the differences in estimated iron core spin rates and an- 
gular momentum profiles in the literature and keeping in mind 
that the study of stellar evolution with rotation is still in its in- 
fancy, one should be cautious in accepting any of the recent 
estimates as a final answer Too many open issues remain and 
can only be resolved by truly multi-dimensional stellar evolu- 
tion studies in the future. Until then, rotating ID evolutionary 
models are the only means we have of estimating the angular 
momentum distribution and spin rate of iron cores at the onset 
of core collapse. During the subsequent collapse evolution, 
the core spins up due to angular momentum conservation and 
— depending upon its initial distribution and amount of angu- 
lar momentum — becomes oblate due to rotational flattening. 
If the rotating core deviates significantly from spherical sym- 
metry, a multi-dimensional treatment becomes indispensable 
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to obtain reliable estimates of neutron star birth periods. 

To date, only a small number of multi-dimensional su- 
pemo va studies have addressed the neutron star spin ques- 
tion. iFrver & Heged ( l2000h performed smooth-particle hy- 
drodynamics simulations with a gray, flux-limited diffusion 
scheme for radiative transfer of the ne utrinos. They used 
the rotating progenitor model E15B of iHeger et al.l ( 120001) 
and estimated an initial spin period of the PNS (at the end 
of their simulation and assuming rigid rotation) on the or- 
der of 100 ms. It is, however, not clear how they defined 
the extent of the PNS. Taking the angular momentum con- 
tained in their PNS and following Hege_Let^aL (2000 ) they ob- 
tained a neutron star spin period of ~2 ms, even for a pro- 
genitor model that included an a-disk-like prescription for 
turbulent-viscous angular momentum redistribution. In their 
recent gr ay 3D s mooth-particle radiation-hydrodynamics cal- 
culation, iFrver & Warren ( 2004) estimated neutron star spin 
periods by assuming that the angular momentum of the inner 
1 M0 is conserved as the PNS cools and con tracts to a neu- 
tron star. For model E15A jHeger et al.l2000h they estimated 
in that way a neutron star period of 0.91 ms and for the slow 
"magn etic" 15 M(n mo del ofiHeger et al. (2003b) a period of 
1 7 ms . lOtt et al.l (l2004h have performed purely hydrody namic 
simulations with a detailed nuclear equation of state and real- 
istic progenitor models and focused on the impact of rotation 
on the gravitational wave signatures of stellar core collapse. 
They argue that if gravitational waves of a core-collapse event 
were observed with sufficient precision it would be possible 
to infer details about the rotational configuration from the ob- 
servations. 

Using flie code VULCAN/2D, IWalder et all (l2005h per- 
formed a series of 2D multi-group, multi-species flux- 
limited diffusion radiation-hydrodynamics supernova sim- 
ulations with the 1 IMq (ZAMS) progenitor model of 
IWooslev & WeaveJ (Il995h . which they artificially forced to 
rotate according to a rotation law that enforces constant an- 
gular velocity on cylindrical shells. Their focus was on the 
anisotropics in the neutrino-radiation field and in neutrino 
heating due to the rotational flattening of the postbounce core 
and they did not systematically analyze the final rotational 
configu ration of their PNSs. In their ID supernova simu- 
lations, [Thompson et al.l (l2005h included the effects of rota- 
tion approximately and they studied the action of viscous pro- 
cesses in dissipating the strong rotational shear profile pro- 
duced by core collapse in a range of progenitors and for dif- 
ferent initial iron core periods. They showed that for rapidly 
rotating cores with postbounce periods of < 4 ms, viscosity 
(presumably by magnetic torques, e.g. the magneto-rotational 
instability [MRI] or magneto-convection) can spin down the 
rapidly rotating PNSs by a factor of ^ 2 - 3 in the early post- 
bounce epoch. 

This paper is the first in a series of papers in which we plan 
to make progressively better estimates of neutron star birth pe- 
riods and rotational configurations, connecting supernova the- 
ory with pulsar science. Here, we present the first systematic 
investigation of the mapping between initial iron core spin and 
the rotational configuration of neutron stars at birth carried 
out with state-of-the-art one- and two-dimensional radiation- 
hydrodynamics codes. We include in our model suite the 
nonrotating lIMfT^, ISMr?. and 2OM0 progenitor models of 
IWooslev & Weave a (fT995l) the r otating 15Mq and 2OM0 
progenitors of iHegeretaD ( I200CI) and the recent ISM© ro- 



tating progenitor of iHeger et al.l ( l2005h that has been evolved 
with angular momentum redistribution by magnetic fields. We 
systematically investigate the dependence of the PNS spin 
on progenitor mass and structure and compare the PNS ro- 
tational configurations resulting from the nonrotating progen- 
itors which we force to rotate according to a prescribed rota- 
tion law to those resulting from models that have been evolved 
with a one-dimensional treatment of rotation. 

In ^ we describe the supernova codes SESAME and 
VULCAN/2D which we employ in this study. Section pre- 
views the progenitor model suite on which we are relyingand 
discusses the rotation laws and mappings used. Section@]de- 
scribes our calculations in one and two spatial dimensions. In 
^ we compare one-dimensional and two-dimensional simu- 
lations and assess the quality of the rotating one-dimensional 
models. In |6l based on the PNS structures that we find in our 
simulations, we estimate the final cold neutron star spin and 
discuss spindown mechanisms. Finally, in ^we summarize 
our results and discuss them critically. 

2. COMPUTATIONAL METHODS 

2.1. 1-D, Spherically Symmetric Simulations 

Our spherically symmetric models are computed with the 
radiation -hydrodynamic s algorithm SESAME^, desc r ibed in 
detail in Burrows et al. (2000) and Th ompson et a l. (2003j). 
The hydrodynamics scheme is Lagrange an, Newtonian, and 
explicit. It is coupled to the EOS of iLattimer & Swestvl 
(TiggT]) (the t.seos), using the implementation described in 
[Thompson etal] (l2003h . The comoving Boltzmann equation 
for neutrinos is solved imphcitly and to order v/c (where v is 
the matter velocity). 

[Thompson et al.l(l2005h have extended SESAME to approx- 
imately include the effects of rotation, as well as angular mo- 
mentum transport by viscosity and local viscous heating. The 
strong negative shear profile (dVL/ dlnr) generated by core 
collapse implies that much of the region interior to the su- 
pernova shock is unstable by the generalized Solberg-H0iland 
criterion and, thus, to the magneto-rotational instability (MRI; 
Balbus & Hawley 1991, 1994; see also Akiyama et al. 2003). 
Under the assumption that magnetic torques provide the dom- 
inant viscosity in the PNS, by equating the shear and Maxwell 
stresses one derives that the shear viscosity should scale as 
^ ~ h^Q,, where h is the pressure scale height. Following 
Thompson et al. (2005), we introduce an a-parameter (in 
analogy with studies of accretion) and take ^ = ah^ft. With 
this form for the viscosity the viscous heating rate scales as 
n^. The viscous timescale is Tvisc ^ {afl)~\r / h)^ . If the 
MRI does operate in PNSs at birth, any seed magnetic field 
grows exponentially on a timescale ~ il"'. Although the sat- 
uration magnetic field strength is somewhat uncertain, to first 
approximation the toroidal magnetic field grows until Bsat, ^ 
(47rp)i/Vfi, which is of order - lO"" G for p = 10'^ g cm-^ 
r ~ 40 km, and a protoneutron star spin period of ~ 10 ms 
just after collapse (Akiyama et al. 2003). The poloidal mag- 
netic energy density saturates at a value approximately an or- 
der of magnitude smaller than B^^^ ^/Stt (e.g., Balbus & Haw- 
ley 1998). 

^ Spherical Explicit/Implicit Supernova Algorithm for Multi- 
Group/Multi-Angle Explosion Simulations. 



4 



All rotating ID models (:1D suffix in the model name) are 
evolved with 500 mass zones, encompassing the progenitor 
iron core. Twenty energy groups are employed in the radia- 
tion transport per neutrino species. Models that include vis- 
cosity have an additional "v" added to their model name (see 
Table All viscous models assume a = 0.1. Viscosity is 
turned on at bounce and is assumed to act at all radii interior 
to the shock radius. 

2.2. 2-D, Axisymmetric Simulations 

Our axisymmetric core-collapse supernova models are 
evolved using the 2D Newtonian, multi-species, multi-group, 
flux-limited (MGFLD) variant of the multi-species, multi- 
group, multi-angle, time-dependent radiation-hydrodynamics 
code VULCAN/2D (Livne et al. 2004). This va riant has pre- 
viousl y been used and described in the study of Walder et alJ 
( 120051) . Its hydrodynamics module (lLivne.l993«) uses an Ar- 
bitrary Lagrangean-Eulerian scheme with remap and a scalar 
von Neumann-Richtmyer artificial viscosity scheme for shock 
handling. W e employ the LSEOS, u sing the implementation 
described bv lThompson et alJ J2003h . 

Since VULCAN/2D uses an Eulerian grid, the specific an- 
gular momentum is advected with the flow in the same man- 
ner as linear momemtum components. In so doing, we main- 
tain global angular momentum conservation by construction. 
Note that the axis in cylindrical coordinates is a singularity 
and, as such, is prone to slightly larger errors than can be ex- 
pected elsewhere on the grid. However, the actual volume of 
the cells nearest the axis is very sm all and the errors do not 
affect the overall flow. In the past dWalder et al.ll200.^ we 
have estimated such departures near the singularity for rotat- 
ing models to be no more than ^--^10% in any flow variable, 
and to be much smaller elsewhere. 

Using the MGFLD variant of VULCAN/2D allows us to 
perform an extensive study that encompasses the evolution 
of many rotating models to times up to ~300 ms after core 
bounce. However, one should bear in mind that MGFLD is 
only an approximation to full Boltzmann transport and dif- 
ferences with the more exact treatment will emerge in the 
neutrino semi-transparent and transparent regimes above the 
PNS surface. Inside the neutrinosphere of the PNS the two- 
dimensional MGFLD approach provides a very reasonable 
representation of the multi-species, multi-group neutrino ra- 
diation fields. T he flux l i miter used in this stu dy is a vector 
versio n of that of lBruennI jl985h . as discussed in lWalder et alJ 
j2005h . In the current version of VULCAN/2D we leave 
out velocity terms, such as Doppler shifts, but include ad- 
vection in the transport sector. Such velocity terms are 
deemed to be of relative importance by the Garching group 
who find neutrino-driven explosio ns whe n turning them off 
(Hampp & Janka 2002; Buras et al1 l2003i 20051. Note, how- 
ever, that velocity terms in Eulerian transport are different 
from the corresponding terms in the comoving frame and that 
statements about their relative importance are very frame- 
dependent. We also do not include energy redistribution by 
neutrino-electron scattering that may affect the size of the ho- 
mologous core at bounce by ^10%, but is otherwise quite 
subdominant. 

A key feature of VULCAN/2D is its abiUty to deal with ar- 
bitrarily shaped grids while using cylindrical coordinates in 
axisymmetry. The grid s etup we u s e for this study is simi- 
lar to the one depicted in lOtt et all J2004h . All our models 



are evolved on the full 1 80 degrees of t he symmetry domain . 
Those models from the simulations of IWalder et alJ i2005t) 
are evolved with 80 evenly spaced angular zones and 130 
logarithmically dist ributed radial zones out to 2000 km. All 
IWalder et alJ J2005ft models are evolved with 8 neutrino en- 
ergy groups for each neutrino species. Our previously unpub- 
lished models are evolved without exception with 16 energy 
groups per neutrino species, 120 angular and 220 radial zones 
out to 3000 km, and, for very massive presupernova cores 
(namely models s20, E15A and E20A; see out to 5000 
km. We have performed a resolution comparison simulation 
with 2/3 the angular and radial zones for a rapidly rotating 
model and find that, despite local qualitative and quantitative 
differences in the flow, integral and global observables, such 
as the moment-of-inertia weighted mean period (see ^ re- 
main practically unchanged, at least at times up to ~150 - 
200 ms after bounce. At the resolution we use in this study, 
VULCAN/2D conserves energy to better than ^0.7% near 
bounce and ~0.4% on average in terms of AE/Egrav At the 
end of the simulation AE/Eg,-av is ^^0.2%. 

3. INITIAL MODELS 

We use two sets of realistic supernova progenitor models 
from stellar evolutionary calculations as initial data for our 
simulations. The first set encompasses models sll, sl5, and 
s20 (ZA MS mass of 1 1 Mf;^ 15 Mp , 20 M©; solar metal- 
licity) of IWooslev & Weaven ( Il995h . These models were 
evolved in spherical symmetry until the onset of iron core col- 
lapse. In Fig. |2j we compare the matter density profiles of 
the presupernova models. Note that models sll and sl5 have 
nearly identical structures out to ^1000 km. The s20 model 
has a shallower density profile and a lower central density, but 
a more massive iron core (~1.9Mq, see Table^. 

On mapping to our simulation grids we introduce rotation 
using the rotation law. 



2n -1 



(1) 



where 51o = 27rPQ' is the angular velocity corresponding to 
the initial central rotation period Pq. A is a free parameter 
that determines the degree of differential rotation in our mod- 
els. For r < A, the initial rotation profile is roughly solid- 
body. According to the Poincare-Wavre theorem, degenerate 
rotating fluid bodies have constant speci fic angular momen- 
tum on cylinders (see, e.g., lTassoull978i) . Hence, we choose 
r in eq. Q to be the cylindrical radial coordinate ("cylindri- 
cal rotation") of the two-dimensional models that we evolve 
with VULCAN/2D. In spherical symmetry, the approximate 
treatment of rotation goes along with the assumption of con- 
stant angular momentum on spherical shells ("shellular ro- 
tation") and we interpret r as the spherical radial coordinate 
in all models evolved in spherical symmetry with SESAME. 
To compare "shellular" and "cylindrical" rotation laws, one 
of our two-dimensional models is set up to rotate initially 
with constant angular momentum on spherical shells. Fig- 
ure^demonstrates that the above rotation law reasonably ap- 
proximates therotationaLgrofiles of the rotating presupernova 
models of iHeger et alJ ( I2OO0I) if one chooses A^IOOO km. 
Except for two models in which we investigate the effect 
of A ^ 00, all other models of this set have A=1000 km. 
We name our models according to the following conven- 
tion: [progenitor name]A[in km]P[Po in s][:evolution type]. 
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For ex ample, sllA1000P8.00:lDv is a IWooslev & Weaved 
J 19951) 1 1 Mq model with A = 1000 km and Po = 8.00 s that 
was evolved in one-dimension with viscous dissipation turned 
on after core bounce. Other evolution types are :1D (one- 
dimensional, no viscous dissipation), :2D (two-dimensional, 
"cylindrical rotation"), 2Ds (two-dimensional, "shellular ro- 
tation") and 2Dh (two-dimensional, "cylindrical rotation," 
purely hydrodynamic, no neutrino transfer). Table [2 details 
our choices of progenitor model, value of A, Pq, and evo- 
lution type for this model set. In particular, we choose the 
range of Pq so that Keplerian rotational velocities are not 
attained anywhere during the numerical evolution, that is 
f2(r,r) < r^Kepiei = \JGM /r^ for all r and t. In terms of the 
widely used rotation parameter /?, which is defined as the ratio 
of rotational kinetic to gravitational potential energy (T/|W|), 
for this series the precollapse values range from 7 x 10"''% up 
to 0.88% for the fastest model. 

The second set of progenitors ("Heger models" in the fol- 
lowing) whi ch we consider consists of mode ls E15 A, E20A 
ileger et ai1 l2000h and ml5b6 CHeger et al1 l2005t) . These 
progenitor models were evolved to the onset of iron core col- 
lapse with the approximate treatment of rotation described by 
Heger et al. (2000), but centrifugal forces were included only 
until the end of core carbon burning. Model E15A and model 
E20A are relatively fast rotators with initial central periods of 
^1.5 and 2 seconds and values of (3 = T/|W| of ^0.5% and 
0.3%, respectively. In model ml5b6, angular momentum re- 
distribution by magnetic processes was included. This led to 
slow and more rigid rotation of its iron core, as can be inferred 
from Fig.[n In Fig.|2l we compare the density profiles of the 
"Heger models" with those of the nonrotating sll, sl5, and 
s20 models. The slow ml5b6 has very similar structure to the 
nonrotating sll and sl5 models. In the 15M0 E15A model, 
rotational effects lead to a considerably more extended iron 
core with lower central density and shallower density strati- 
fication than in the nonrotating sl5. Model E20A's structure, 
however, does not deviate significantly from that of the non- 
rotating s20. 

The models are mapped onto our computational grids ac- 
cording to the same prescription presented for the first set: 
"shellular rotation" in ID, "cylindrical rotation" for most two- 
dimensional models. In naming the models belonging to 
this set, we follow the conventions mapped out for the first 
set, but leave out the parts concerning paramete rs A and P: 
E15A:2Ds, for example, is progenitor E15A of iHeger et al.l 
(|2000) evolved in axisymmetry, but with "shellular rotation." 
Table[niists the initial model parameters for E15A, E20A, and 
ml5b6. 



4. RESULTS 

The focus of this study is on the periods and rotational pro- 
files of the PNSs as obtained from our one-dimensional and 
two-dimensional radiation-hydrodynamics simulations. For 
more in-depth discu ssions of the supernova physics, we refer 
the reader to Thomp son et al.l fcOOS) fo r the one-dime nsional 
simulations with SESAME and to WakLer et al.l (120051) for the 
two-dimensional MGFLD VULCAN/2D simulations. 

In the discussion of our one-dimensional simulations we 
limit ourselves for conciseness to tho se performed with the 
sll model (11 Mq at ZAMS) of iWooslev & Weaved (IT995h 
while we describe and compare in detail simulations with a 



wide variety of progenitor star models in the discussion of 
our two-dimensional calculations. 

4.1. One- dimensional Simulations 

An interesting consequence of the self-similar, quasi- 
homologous nature of the inner core's collapse is the preserva- 
tion of its initial rotational profile in the PNS from the precol- 
lapse state. In other words, if the inner core is initially rotating 
rigidl y (which is likely to be the case; see lHeger et al.ll2000l 
l2005l andFig.m. then the inner PNS will be in solid-body ro- 
tation as well. There is, however, a region in any PNS that is 
rotating differentially — independent of the angular momen- 
tum distribution in the progenitor core. Due to angular mo- 
mentum conservation the core spins up during collapse. Since 
the inner parts of the precollapse iron core collapse to smaller 
radii than the supersonically collapsing outer parts, differen- 
tial rotation is a natural consequence of collapse. However, 
solid-body rotation is the lowest-energy state of a rotating 
body and viscous processes can drive the differentially ro- 
tating PNS into rigid rotation as it evolves, cools and con- 
denses to its final configuration. As iThompson et alJ (l2005h 
have pointed out, early and efficient postbounce dissipation 
of the shear energy stored in the differential rotation of the 
PNS by viscous processes (e.g., the MRI) can significantly 
affect the dynamics and might result in a successful explo- 
sion. To estimate the rotation periods our PNSs would have 
if they were in solid-body rotation, we introduce a moment of 
inertia-weighted mean angular velocity: 



X^oCMi-)^(r)ridV 

/*oCp(r)ridV 



(2) 



where r±_ is the distance from the rotation axis at any given 
r, r = |r| is the spherical radius, R is the outer radius of the 
computational domain and p{f) is the matter density. In the 
one-dimensional case, we use the spherical radius for r^. We 
now define two mean rotation periods. The definition of the 
function C depends on the individual mean being considered. 
The first mean period we define is 



P^i2 = 27r(17(C))-' , where 
r)=/l ifp(r)> lO'^gcm- 



C(r) = 



otherwise 



(3) 



Ppi2 is the period that all material with p > lO'^ g cm"'^ would 
have if it were in solid-body rotation. We use this quantity 
to describe the postbounce spin evolution of the PNS. It is 
particularly useful because it captures the effects of accretion 
and contraction to higher densities and is meaningful in both 
one and two spatial dimensions. We adopt the density cut 
of 10'^ g cm"^, roughly marking the characteristic density at 
which neutrinos decouple from the fluid, to define the "edge" 
of the PNS for analysis purposes. Clearly, a PNS does not 
have a sharp boundary. However, most of its material has 
densities above the assumed density cut. The second mean 
period is 



PM=27r(17(C))"' , where 

f 1 ifM(r) <M 
^ ' 10 otherwise 



(4) 



Pm is the period that the material interior to a given mass 
coordinate M would have if it were in solid-body rotation. 
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An unambiguous definition of this quantity is only possible in 
one spatial dimension, because, in two dimensions one cannot 
assume that the matter distribution remains spherically sym- 
metric — rotation and postbounce convection break spherical 
symmetry in any two-dimensional model. 

The evolution of the postbounce mean period Ppi2 as de- 
fined above is shown in Fig.|3lfor a variety of models of the 
si 1 model series. Shown in solid lines is the mean period for 
models without active viscous dissipation (:1D models). In 
dotted lines we plot the corresponding models in which we 
include viscous dissipation (:lDv models). Let us first fo- 
cus on the evolution of Ppi2 for models without the inclusion 
of viscous effects. During collapse, P^n remains undefined 
until the density reaches the threshold value of 10'^ g cm"^^. 
Initially it increases as more slowly rotating material achieves 
a density of 10'~ g cm"-' and reaches its maximum just ^1 ms 
before core bounce, when the lion's share of the inner core has 
crossed the density threshold, but has still not spun up to its 
bounce angular velocity. At bounce, Ppi2 experiences a local 
minimum and then increases again as the core re-expands to 
a quasi-equiUbrium configuration. 

Model si 1A1000P1.25:1D (shown in black in Fig.|3} is our 
fastest rotator, with an initial central period of 1.25 s and (3 = 
T/\W\ ^0.9%. It undergoes a core bounce altered by centrifu- 
gal forces and forms a quickly spinning PNS with a central 
period of ^^1 ms. As displayed in Fig.|3] its mean period Ppi2 
stays almost constant at ~3 ms until the end of the calculation. 
Because of strong centrifugal support, its protoneutron star is 
large and, even though matter is accreting, is only slowly con- 
tracting. Furthermore, model sllA1000P1.25:lD is the only 
one-dimensional model without viscous dissipation that ex- 
hibits an explosion — though a weak one — and less matter 
is accreted onto the PNS than in models that do not explode. 
Model sllA1000P2.00:lD is the next slowest rotator in our 
sll model series. Its core bounce is dominated by nuclear 
repulsive forces, but centrifugal forces still play an important 
role. Its Ppi2 increases slightly after bounce, reaches a local 
maximum at a time of ^^30 ms after bounce, and decreases 
from thereon to a finaF"* Ppi2 of 3.5 ms. 

With increasing initial central rotation period, the influ- 
ence of centrifugal forces becomes less and less relevant. As 
can be inferred from Table |2] the initial value of [3 is al- 
ready as low as 0.25% for model sllA1000P3.00. In Fig. 13 
all models with initial central period above 3 seconds show 
the same qualitative behavior: a short increase of Ppi2 af- 
ter bounce and then a monotonic decrease as more mate- 
rial is crossing the 10'~ g cm"^ density threshold and piling 
up on the PNS. The slowest model evolved without viscous 
dissipation, sllA1000P10.47:lD, has a final (~500 ms after 
bounce) mean period of ^14 ms. Model sllAooP10.47:lD 
(magenta graph in Fig.|3ji has a rotation parameter A (eq. Q) 
set to infinity, forcing the initial model to rotate like a 
solid body. Compared with model sllA1000P10.47:lD, 
si 1 AooPlO.47: ID has significantly more angular momentum 
in the outer regions of the core. This difference has a major 
impact on the PNS's rotation and leads to a ~50% smaller 
final mean period. 

Dissip ation of the shear energ y stored in differential rota- 
tion is, as lThompson et al]i2005t) report, most efficient in fast 



rotating models. This is supported by the evolution of our sll 
models that include the effects of viscous energy dissipation 
and angular momentum redistribution after bounce. In Fig.|3| 
we plot in dotted lines the evolution of Ppi2 for our : IDv mod- 
els. For model si 1A1000P2.00, the :lDv variant begins to 
deviate significantly from the model without viscous dissi- 
pation at about 40 ms after bounce. While the :1D model's 
mean period slightly decreases, Ppn of the :lDv model in- 
creases strongly as long as the angular momentum redistribu- 
tion mechanism provided by viscosity remains efficient, that 
is, as long as there is considerable shear energy in differen- 
tial rotation that can be dissipated and as long as the vis- 
cous timescale is short with respect to the Kelvin-Helmholtz 
timescale. Eventually, at ~300 ms after bounce the entire 
PNS is in solid-body rotation and Ppi2 remains constant, or 
decreases slightly, as the PNS contracts. All fast models for 
which we have performed :lDv evolution (si 1A1000P2.00, 
S11A1000P2.34, S11A1000P3.00, sllA1000P4.00) experi- 
ence succe ssful explosions (for det ails on the explosion mech- 
anism, see iThompson et ai1l2005h . Since rotational energy 
scales with fi", there is, for the same qualitative rotational 
profile, much less energy stored in the differential rotation 
of slow models than in fast models. This is why slow one- 
dimensional rotators (Pq > 5 s) do not explode even with 
the additional heating provided by viscous dissipation. Fur- 
thermore, this comports with the fact that the change in Ppi2 
brought about by the inclusion of viscous effects operating 
on differential rotation becomes smaller with increasing ini- 
tial period. (Recall tha t the viscous timescale scales with O"' ; 
IThompson et alJ2005[ ') 

The global ratio of rotational kinetic energy T to gravita- 
tional potential energy |W|, the rotation parameter /3 = T/|W|, 
is another useful and important quantity whose evolution we 
follow for all our models. If the value of [3 passes some 
critical threshold value, nonaxisymmetric rotational fluid in- 
stabilities may set in and lead to large-scale deformation 
of the core. In the approximation of MacLaurin spheroids 
(uniform density, incompressible, rigidly rotating equilibrium 
configurations), triaxial instabilities may grow if /3 > 14% 
or P > 27% for sec ular and dynamical inst abilities, respec- 
tively~iTassoul"197^. Given recent results by Ott et al.'('200S 
(see a lso Centrella et al. 2001; Shibata et al. 2002; Saijo et a2 
1200 3'; Shibata et al. 2003), it is likely that the above canonical 
threshold values for f3 do not hold for differentiall y rotating 
postbo unce cores with a realistic equation of state. lOtt et alJ 
J2005h have found that a dynamical m = 1 -dominated spiral 
instability may set in at a /3 as low as 8%. Nonaxisymmetric 
instabilities are of particular interest because they may lead 
to the emission of strong gravitational waves and to angular 
momentum redistribution and spindown of the central core. 

Figurel^shows the evolution of (3 for selected : ID and : IDv 
models of the sll model series. sllA1000P1.25:lD is the 
fastest rotator in this series. Initially, [3 is 0.88% for this model 
and it reaches a local maximum of 14.6% at bounce. After 
bounce, the relaxing core, abetted by centrifugal forces, over- 
shoots its equilibrium configuration, leading to a pronounced 
local minimum in (3 of 12% about 3 ms after bounce Im- 
mediately following the initial postbounce relaxation, f3 be- 
gins to increase and does so throughout the rest of the evolu- 
tion, reaching a /?/ (the / indicating "final") of 18.4% at the 



^ In the following — unless otherwise stated — "final" means "at the end 
of the numerical evolution." 



Due to the range in time we are showing in Fig.|4] this local minimum is 
not clearly visible. 
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end of our numerical evolution. Given the approximate con- 
stancyof the mean period P^n of model si 1A1000P1.25:1D 
(Fig. 13 the above might seem surprising. However, /3 is a 
global quantity; the integrals for T and W encompass the en- 
tire grid and not just material above a density cut. Hence, 
the linear (in time) increase in global /3 indicates that the pro- 
toneutron star is slowly beginning to contract and spin up, but 
also that more and more material from larger radii is falling 
in, spinning up and settling onto the PNS surface. Both T and 
|W| increase, with T dominating, because of its higher-power 
dependence on the radial coordinate. 

The f3 evolution for the rest of the sll models without vis- 
cous dissipation is essentially the same as the one just de- 
scribed for model sll AlOOOPl .25: ID. With decreasing initial 
f3, (3 at bounce and /3/ decrease. Model sllAooP10.47:lD, 
whose core is initially in solid-body rotation, stands out. 
Even though its initial period Pq is the same as model 
sllAlOOOPlO.47: ID'S, its initial value of (3 is more than four 
times larger (see Table |3 than sllA1000P10.47:lD's. As 
shown in Fig.0] this is reflected in its /3 evolution up to bounce 
and, most apparently, in the steeper increase in the postbounce 
phase. 

Two of our models from the sll model series, namely 
sllA1000P1.25:lD and si 1A1000P2.00:1D surpass /3 = 
8% during their numerical evolution. They also exhibit 
strong differential rotation (Fig. |5} and mig ht experience th e 
m=l-dominated one-armed spiral instability (lOtt et al.l2005l) . 
sllA1000P1.25:lD even reaches (3 > 14% and may become 
unstable towards a "classical" bar-like instability on a secular 
timescale. 

The evolution of /3 for models with viscous angular mo- 
mentum redistribution is quite different (Fig.|4j. A s soo n as 
the viscous terms are switched on after bounce ( ^2.1> . for 
:lDv models with Pq < 4 s (red, green and brown dotted 
lines in Fig. |^ /? begins to decrease. Shear energy stored 
in the differential rotation of the PNS is dissipated away, 
leading to solid-body rotation at lower angular velocities 
an d, hence, lower rotation al kinetic energies. As discussed 
bv lThompson et alJ J2Q05I) . models si 1 A1000P2.00: IDv and 
si lA1000P3.00:lDv experience energetic explosions that re- 
move a significant fraction of the stellar envelope. There- 
fore the PNS masses remain somewhat smaller than for 
sllA1000P4.00;lDv, which explodes less energetically and 
at later times, allowing more matter to accrete through the 
stalled supernova shock and onto the PNS (Fig.|5). In model 
sllA1000P4.00:lDv, /3 decreases for the first 100 ms and 
then increases again. This is, in part, because of the ^ il^ 
dependence of the viscous dissipation and, in part, because of 
its weak and later explosion that results in more mass being 
accreted than in the two faster models discussed. For even 
slower models, the evolution of /3 is similar to those for the 
cases without viscosity. Model si lAooP10.47:lDv departs 
from the :1D evolution at late times when high-angular mo- 
mentum material reaches small radii and the energy stored in 
differential rotation becomes large. Viscous dissipation of that 
energy then leads to a smaller final /3 than in the :1D model 
without viscous dissipation. 

In Fig. |5] we present the final period (Fig. (5^) and an- 
gular velocity profiles (Fig. of a representative subset 
of the sll model series calculated with SESAME. Profiles 
of models with viscous dissipation and angular momentum 
redistribution are drawn as dotted lines, profiles of mod- 



els without viscous effects are drawn as solid lines. Model 
sllA1000P1.25:lD, drawn in black, has the shortest initial 
central period (Pq = 1.25 s) which maps to the shortest fi- 
nal central period of Pq,/ = 1.56 ms. As discussed ear- 
lier in this section, the PNSs are in solid-body rotation in- 
side 8-10 km (or ^0.5 - 0.8 Mq in mass coordinate). 
At greater radii, the PNSs quickly become differentially ro- 
tating. For model si 1A1000P1.25:1D, the angular velocity 
drops from --2500 rad s"' at 20 km to --900 rad s"' at 50 km 
and ^^330 rad s~' at 100 km, corresponding to a Ail/Ar of 
53 rad s"' km"' from 20 to 50 km and 1 1 rad s"' km"' from 
50 to 100 km (be aware of the log-log scale in Fig.|5}- Again, 
one has to keep in mind that model sll AlOOOPl. 25: ID ex- 
plodes. The explosion removes material that would otherwise 
accrete onto the PNS and pile up on the quasi-incompressible 
inner region above radii of 15 - 30 km. This is reflected 
in the features visible in the graphs of :1D models with 
Po > 2 s in the range of radii mentioned. Models in which 
we include viscous dissipation and angular momentum trans- 
port exhibit solid-body rotation out to much larger radii. 
Model sllA1000P2.00:lDv, for example, rotates rigidly out 
to ^23 km or 0.9 Mq in mass coordinate while model 
sllA1000P2.00:lD exhibits solid-body rotation out to only 
~10 km or 0.42 Mq^. Viscous dissipation of shear energy 
stored in differential rotation also leads to a spindown of the 
PNS. sllA1000P2.00:lD's final central period is 1.95 ms 
whereas si 1 A1000P2.00: IDv spins with a final central period 
of 6.7 ms, which is more than 3.5 times slower Due to the 
dependence of rotational energy and the 51^ dependence 
of the dissipation rate, the relative period increase brought 
about by viscous dissipation/angular momentum transport de- 
creases quickly with increasing Pq. For si 1 A1000P3.00 (cen- 
tral periods: 2.74 ms and 6.16 ms for the :1D and :lDv 
models, respectively), the ratio is down to about 2.2 and for 
sllAlOOOPlO.47 (central periods of 9. 19 ms and 16.0 ms for 
models :1D and :lDv, respectively), it is approximately 1.7. 
Of course, the just quoted ratios depend on the simulation 
time after bounce, which determines how much time we have 
allowed the dissipative process to spin down the core and how 
much material and angular momentum is accreted during the 
simulation. Model si lAooP 10.47 stands out. It contains sig- 
nificantly more angular momentum than its counterpart with 
A = 1000 km. Therefore, it spins up to shorter central periods 
during collapse and during the postbounce evolution. 

For key models of the sll model series, the final moment- 
of-inertia-weighted mean period Pm , as defined by eq. 0, 
is shown in Fig. |6^. Note that Pm gives the mean period of 
material inside a mass coordinate M, assuming solid-body ro- 
tation. Hence, the graphs should not be interpreted as profiles 
in the usual sense. The mass coordinate at which the slope 
of each individual model's Pm sharply steepens is a good es- 
timate for the compact remnant mass in models that exhibit 
explosions, si lA1000P2.00:lDv (dotted red), for example, 
explodes early and energetically, leaving behind a lower-mass 
PNS than model sllA1000P3.00:lDv (dotted green), which 
explodes at a somewhat later time. For models that don't ex- 
plode, Pm jumps only at ~1.42M0. The consequences of vis- 
cous dissipation and angular momentum redistribution in the 
:lDv models, spindown and solid-body rotation, are perhaps 
even more obvious in Fig.|6t than in Fig.|5] Also, note that, 
since the :lDv models are in solid-body rotation throughout 

^ See Fig.|6|) for the mapping between mass and radial coordinates at the 
end of each model's numerical evolution. 
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most of the protoneutron star at the end of our calculations, 
Pm for them agrees very well with the actual period profiles. 



4.2. 2D Simulations 

We now discuss the results of the axisymmetric 2D VUL- 
CAN simulations using the methods described in ^2.21 
In contrast to the spherically symmetric simulations, two- 
dimensional, azimuthally-symmetric simulations are able to 
follow the non-spherical effects of rotation in core collapse. 
The degree of oblate deformation can be described by the ra- 
tio of the PNS's polar and equatorial radii on isodensity sur- 
faces, the so-called axis ratio. FigureQdisplays snapshots of 
the density at ~'200 ms after bounce for a set of representa- 
tive 2D models. Models sllA1000P2.34 and E15A are the 
fastest rotators, having axis ratios below 1:2. Clearly, a two- 
dimensional treatment is a must for capturing the full dynam- 
ics of these rapid rotators. With increasing initial iron core 
period (Pq) the PNS axis ratio increases towards one. ml5b6 
and the two variants of si 1A1000P10.47:2D exhibit almost 
no rotational flattening at all and one might expect a one- 
dimensional treatment to yield good results. But it is not only 
rotation that is modeled much more realistically in real 2D 
simulations. Neutrino-driven convection naturally develops 
in the region behind the stalled supernova shock and leads to 
the interesting substructure seen in Fig.0 The purely hydro- 
dynamic model, sll A1000P10.47:2Dh, experiences a prompt 
explosion and Fig. does not show the imprint of the vig- 
orous convection present in the models that are evolved with 
adequate neutrino treatment and radiative transfer Note that 
the small dents and spikes visible along the rotation axis in 
most of the panels are related to the imperfect symmetry axis 
treatment in the version of VULCAN/2D with which most of 
our models are evolved. The purely hydrodynamic model has 
been evolved with an updated version of the code and exhibits 
minimal axis artefacts. 

Because of the computational difficulties that accompany 
long-term, two-dimensional radiation-hydrodynamics evolu- 
tions, we have run our two-dimensional models to postbounce 
times of only 200 - 300 ms and display profiles and two- 
dimensional snapshots at 200 ms after bounce. At that epoch, 
the PNS has already reached a compact state and the largest 
fraction of its final mass has already settled inside a (spher- 
ical) radius of ~50 km. Figure |8] shows the mass contained 
inside a given spherical radius for selected models at 200 ms 
after core bounce. The graphs certainly ignore the two- 
dimensional nature of the calculation, but give a good (and 
for the slow models, almost accurate) handle on the mass- 
radius relationship. The dashed lines in Fig.|S]represent mass 
profiles of models that are evolved adiabatically and exhibit 
prompt explosions, blasting away a large fraction of the outer 
iron core and the stellar envelope. Hence, their PNS masses 
are small compared with those of the MGFLD models that 
might (or might not) explode at later times after much of the 
outer core has accreted onto the PNS. Note that slower mod- 
els have more compact PNSs than faster models. This is, of 
course, expected since rotational support is greater in faster 
models and leads to more expanded postbounce configura- 
tions. 

In Fig. |5] we depict the equatorial angular velocity and 
period profiles (panels (a) and (b), respectively) at 200 ms 
after bounce for the two-dimensional sllAlOOO model 
series. For comparison we also display graphs of model 
sllAooP10.47:2D (magenta; initially in solid-body rotation) 
and of two purely hydrodynamic models (dashed graphs). 
Note that the graphs start at a radius of 3 km to excise axis 
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artefacts. The aim of Fig. |9] is to display the systematic 
dependence of the protoneutron star rotational profile on the 
initial iron core spin rate. For this, we focus in Fig. |9]on 
the angular velocity and period profiles of the MGFLD 
models sllA1000P157.10;2D, si 1A1000P41.69:2D, 
sllA1000P10.47:2D, sllA1000P4.69:2D, and 

si 1A1000P2.34:2D. These models form a sequence in 
initial rotation period from 157.10 s down to 2.34 s. At 
200 ms after bounce, sllA1000P157.10:2D has reached a 
central period of ^^50 ms. The spin period increases by about 
a factor of two from 3 km to 5 km and then stays roughly 
constant, but shows significant substructure, out to 30 km. 
From 30 km to 100 km, the period increases by a factor of 
ten, corresponding to a "period slope" of ~30 ms km"'. We 
call the reader's attention to the region of quasi solid-body 
rotation between roughly 5 and 30 km also present in all 
spherically symmetric models (see ^4.11 and Fig. |5l). This 
region is a direct consequence of self-similar collapse of 
the inner core. The equatorial angular velocity and period 
profiles of the two-dimensional models are, however, much 
less smooth than the profiles of the one-dimensional models 
. This is attributed to the much more complicated, two- 
dimensional dynamics that involve convection and inviscid 
angular momentum redistribution along the rotation axis due 
to convection (Fig. I10> . The hump in the angular velocity 
profile (corresponding to a notable dent in the period profile) 
at about 20 km is a generic feature in the entire sllAlOOO 
series, but is most pronounced in the slow models. It is 
caused by accreting relatively high-specific angular momen- 
tum material that is piling up on the nuclear-density region 
of the PNS. The faster a model is, the more extended its PNS 
(see Fig.|8} and the shallower its density profile. 

The qualitative shape of the angular velocity and period 
profiles of si 1A1000P157.10:2D is mirrored, but of course 
shifted, in the fis (and periods) in faster models (Fig. |9}. 
However, the profiles of faster models are smoother. Con- 
vection in fast models is largely suppressed by rotation in the 
equatorial regions ^ (Frver fc Hegel l20()(it iFrver & WarrenI 
12004: WaMgi"_et_a l- 2005 ) and, hence, does not create as 
much substructure in the equatorial profiles as in slower 
models. In Fig. |9l for models si 1A1000P2.34;2D and 
si 1A1000P10.47:2D, we include equatorial profiles of their 
purely hydrodynamic :2Dh variants. They are plotted with 
dashed lines. Since we do not include any neutrino physics 
and transport in the :2Dh simulations, these models cannot 
contract and spin up to higher after core bounce. They expe- 
rience prompt explosions and fall-back accretion is negligible. 
Keeping this in mind, it is a bit surprising to see how similar 
the hydrodynamic models' profiles are to the profiles of their 
MGFLD counterparts. To conclude our discussion of Fig.|9l 
we call attention to the graphs of si 1 AooP10.47:2D's angular 
velocity and period profile (magenta lines). This model is set 
up in solid-body rotation with an initial spin period of 10.47 s. 
Initial solid-body rotation means that a great part of the iron 
core's angular momentum is located at large radii (the mo- 
ment arm balances the density drop with increasing radius). 
Due to angular momentum conservation during collapse, one 

* Note, however, that for convection to be stabihzed, the rotation rate 
must be large enough that if the calculation had been done in magneto- 
hydrodynamics, the r egion would likely be destabilized by the MRI 
IThompson et aljl2005l> . We therefore expect that when these rapidly ro- 
tating models are calculated in MHD, they will show enhanced (magneto- 
)convection with respect to the hydrodynamical, centrifugally-stabilized cal- 
culations presented here. 



expects the outer core layers to spin up appreciably. This is 
exactly what we see in Fig.|9l The region of high angular ve- 
locity (~600 rad s"') and low period (~10 ms) extends out to 
about 70 km, whereas si lA1000P10.47:2D's angular veloc- 
ity (orange lines) starts dropping (and its period starts increas- 
ing) already at about 20 km. Note that the centermost regions 
of the two models agree relatively well — as they should — 
given the fact that si lA1000P10.47:2D's initial angular ve- 
locity is almost constant throughout its inner core (Fig.Q. 

The fastest model of the series presented in Fig. |9] is 
sllA1000P2.34:2D. Its initial central iron core spin (Pq) of 
2.34 s maps to a central postbounce period (at the 200 ms 
snapshot shown here) of about 1.7 ms, which increases to 
~3.2 ms at 10 km. As shown in Fig.Q its PNS is rotationally 
flattened with an axes ratio below 1:2. In Fig.^| we show the 
angular velocity distribution in si 1 AlOOO models with PqS of 
2.34 s, 4.69 s, 10.47 s, and 41.89 s at 200 ms after bounce. 
Since we are imposing axisymmetry, the angular velocity dis- 
tribution is symmetric with respect to the vertical axis (the 
rotation axis). We superpose flow velocity vectors to convey 
qualitatively the flow dynamics in the PNS. We downsample 
the number of velocity vectors at radii that are greater than 
the approximate shock location of each individual model. For 
the purely hydrodynamic model, the radius beyond which we 
downsample is chosen to be 150 km. The panels show the 
inner 400 km of the computational domain. (For a more de- 
tailed discussion of the multi-dimensional flow dynamics and 
its consequences concerning the sup ernova mechanism, we 
refer the reader to lWalder et alJl2()0a .') All 2D models shown 
rotate in the precollapse stage with constant angular veloc- 
ity on cylindrical shells. As the collapse proceeds, the matter 
falls in almost spherically (more spherically for slowly spin- 
ning models, less spherically for fast spinning models). The 
initial conditions and the collapse dynamics are reflected in 
the elongated (prolate) shape of the angular velocity distribu- 
tion of the PNSs. The fast models (si 1A1000P2.34:2D and 
si 1A1000P4.69:2D) have a notable central spheroidal bulge 
of high angular velocity material, which is formed by high 
specific angular momentum material from the equatorial re- 
gions. 

As in the discussion of our one-dimensional models in 
^4.11 we define a mean, moment-of-inertia-weighted pe- 
riod Ppi2 for all material with density above 10'~ g cm"^* 
(eqs. 121 and |3}. Figure ^3 depicts Ppn as a function 
of time for all two-dimensional sll models. In addi- 
tion to the sll AlOOO MGFLD models, we plot Ppn for 
si lAooP10.47:2D (magenta) and the purely hydrodynamic 
variants si lA1000P2.34:2Dh and si lA1000P10.47:2Dh 
(dashed lines). As one would expect, all models exhibit a 
local minimum of P^n at the time of core bounce. For the 
purely hydrodynamic models, this is the absolute minimum 
in Ppi2 , since they quickly achieve hydrostatic equilibrium 
after bounce. Their periods remain almost constant at post- 
bounce times. Model si lA1000P2.34:2Dh (dashed red) un- 
dergoes core bounce aided by centrifugal forces and shows 
the pronounced postbounce coherent expansion a nd recon- 
tractio n cycles described for these sorts of models in lOtt et all 
(l2()0l . In contrast, si 1A1000P2.34:2D^ its MGFLD vari- 
ant, relaxes from its local minimum in Ppi2 and shows no 
oscillatory behavior at all. We attribute this major difference 
to the following phenomenon. In adiabatic collapse, angu- 
lar momentum conservation during collapse leads to an in- 
crease in the kinetic energy stored in rotation. By its nature 
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rotation acts approximately lik e a gas with F of 5/3 jTassoull 
ll978tlMonchmever et al.ll991l) and in this way effectiyely in- 
creases the pressure support of the collapsing core. At the mo- 
ment of strongest compression at bounce, centrifugal forces 
are maximal and result in the core dramatically overshoot- 
ing its postbounce equilibrium position. Recontraction fol- 
lows, leading to another expansion-recontraction-bounce cy- 
cle. Eventually, enough pulsational energy is dissipated by 
the generation of secondary shocks that after several cycles 
the oscillation is damped. However, including the production 
and radiative transfer of neutrinos, the rebounding core loses a 
large amount of energy in the neutrino burst that occurs when 
the bounce shock wave breaks through the electron neutri- 
nosphere. This, and continuing neutrino loses, lead to en- 
hanced damping of the postbounce expansion-recontraction- 
bounce cycles, far more than experienced in quickly spinning 
purely hydrodynamic models. 

We note that the qualitative shape for a given progenitor 
structure and rotational "scale height" A of the Ppi2 evolu- 
tion shown in Fig. [^is the same for all sllAlOOO mod- 
els. Generally, at any postbounce time Ppi2 is a monotonic 
function of the initial iron core spin. At 200 ms after core 
bounce, the largest Ppi2 is -180 ms (si 1A1000P157.10:2D) 
and the smallest is -4.7 ms (si 1 A1000P2.34:2D). In Table|3l 
we list the mean periods for all models at the end of their 
numerical runs. The magenta graph in Fig. ^2 represents 
si lAooP10.47:2D'sPpi2 evolution. As expected, this initially 
rigidly rotating model's PNS is spinning faster than its coun- 
terpart with A set to 1000 km before collapse. 

As demonstrated for the sll model series in Fig. the 
systematics with the global rotation parameter (5 is similar 
Up to bounce, [3 increases monotonically in all models. The 
faster of the two purely hydrodynamic models (dotted-line 
graphs) reaches the highest {3 at bounce of all models pre- 
sented in Fig.^l Due to its hig her electron fraction Yg and 
higher effective adiabatic F fM artinez-Pinedo et alJ l20()4l) . 
in direct contrast with the more realistic MGFLD variant, 
si 1A1000P2.34:2D, the purely hydrodynamic model has a 
much more massive inner core. Hence, during its "plunge" 
phase just before and at bounce, more material, and in that 
way more angular momentum, reaches small radii. This is 
why purely hydrodynamic models reach higher /3s at bounce 
than their MGFLD diffusion counterparts. While the purely 
hydrodynamic models settle at a practically constant final /?, 
all MGFLD models exhibit a linear postbounce growth of the 
rotation parameter, with a slope that depends on the initial iron 
core spin and on the scale parameter A in the rotation law. 
The increase in [3 is caused by the increase in compactness 
and by accretion of outer-core material. One might expect 
that the slope of j3 will change after a successful explosion 
(not tracked in our models). At the end of their numerical 
evolution none of the models presented in Fig.[2]has reached 
the classical threshold values of B for s ecular or dynamical 
instability. However, as lOtt et al.l (l2005h have demonstrated, 
at P —8% realistic postbounce cores might undergo the low- 
T/|W| instability. This value is likely to be reached by model 
si 1 A1000P2.34:2D at late times. 

Having described the quantitative and qualitative features 
of a neutron star's birth spin and its dependence on the initial 
iron core spin rate for a fixed progenitor structure, we now 
investigate the effect of varying progenitors on the PNS spin. 
We have performed a set of calculations employing the sl5. 



s2 (15 Mf:i and 20 Mp at Z AMS, respectively) progenitors 
of lWooslev & Weave J (i 19951) the E 15A and E20A presuper- 
nova models of iHeger et alJ (l2000h . and the ml5b6 progeni- 
tor of [Heger et alJ (l2005h . The key progenitor characteristics 
of all models are summarized in Tabled In the following 
we present the results of these calculations and compare them 
with those of our sll models. 

In Fig. we contrast the two-dimensional angular ve- 
locity distributions of models with different progenitor iron 
core mass and structure at —200 ms after core bounce. We 
zoom in on the inner 400 km of our computational domain. 
The rotation law parameters Pq and A are prescribed. The 
top right panel of Fig. shows an Vl snapshot of model 
si lA1000P2.34:2Dh, the purely hydrodynamic variant of 
si 1 A1000P2.34:2D. At -200 ms after bounce, this model has 
already undergone a prompt explosion. This explains the out- 
ward pointing fluid velocity vectors in this panel. The region 
of high angular velocity is stretched out along the axis of ro- 
tation (the vertical). This is caused by the prompt explosion, 
which leads to an ejection along the axis of high- 51 material 
from deep within the inner core. 

The lower two panels of Fig. show angular ve- 
locity distributions of models sl5A1000P2.34:2D and 
s20A1000P2.34:2D, respectively. Comparing them to the top- 
left panel which shows model si 1 A1000P2.34:2D, one notes 
that the PNS rotational structure is not strongly dependent on 
progenitor mass; the three models in question have qualita- 
tively and quantitatively very similar angular velocity distri- 
butions at 200 ms after bounce. For the sll and sl5 mod- 
els, this is not surprising, given that their iron core masses 
and density profiles (Fig. |2] and Tabled are very similar to 
each other However, the s20 progenitor model which we 
use in our model s20A1000P2.34:2D has an appreciably dif- 
ferent initial iron core structure. It is more extended, has a 
lower initial central density, and is significantly more mas- 
sive than the sll and sl5 models. Yet, at 200 ms after 
bounce only small differences in the angular velocity dis- 
tribution are apparent. s20A1000P2.34:2D's PNS is more 
massive and more extended. This is reflected in the slightly 
larger central spheroidal bulge of high-angular-velocity ma- 
terial. As can be inferred from the flow patterns in Fig. 
the thermodynamic structure and the resulting fluid dynam- 
ics of models si 1A1000P2.34:2D, sl5A1000P2.34:2D, and 
s20A1000P.23:2D are significantly different. These differ- 
ences are reflected in the angular velocity distribution to only 
a slight degree. 

Figure depicts the angular veloc ity distributions of 
the recent rotati ng "Heger models" jHeger et"al] 120001 
IHe ger et al.ll2005l) at 200 ms after core bounce. We include 
models E15A, E20A, and ml5b6 in our two-dimensional 
model set and for model E15A:2D also perform a simula- 
tion with initial "shellular rotation" (E15A:2Ds; see The 
upper left and upper right panels of Fig. are snapshots 
of the two-dimensional Vl distributions of models E15A:2D 
and E15A;2Ds, respectively. One can easily discern the 
more sphere-like VL distribution of E15A:2Ds at large radii. 
E15A:2Ds' angular velocity in the central region exhibits a 
somewhat more spherical shape than that of its E15A:2D 
counterpart. However, the magnitude of the angular veloci- 
ties observed does not depend strongly on the initial choice 
of rotation law, since the equatorial regions of the precol- 
lapse core that undergo the greatest spin-up during collapse 



11 



have almost identical before collapse (eq.^and the discus- 
sion in 13 ■ The bottom-left panel of Fig. 1141 shows model 
E20A;2D's angular velocity snapshot at 200 ms after bounce. 
Model E20A has an initial core mass and density profile that 
are similar to those of the s20 progenitor (Fig. |2j- Looking 
at Fig. n one notes that E20A's initial ft profile exhibits an 
order-of-magnitudejump at around 1000 km, whereas EISA's 
angular velocity jumps at about 3 times that radius. This dif- 
ferent behavior in the initial il profile and the initially longer 
central period of the E20A progenitor (Table Q would natu- 
rally have led one to anticipate qualitative and quantitative dif- 
ferences in E20A:2D's angular velocity evolution when com- 
pared with that of models E15A:2D and s20A1000P2.34:2D. 
In fact, E20A;2D exhibits an appreciably different angular ve- 
locity profile. At lateral angles between about 25 and 155 
degrees, the angular velocity is almost an order of magni- 
tude lower than in the polar wedges. This can be understood 
as follows: At 200 ms after bounce, the high-angular veloc- 
ity material that was inside a cylindrical radius of ~1000 km 
has in the equatorial regions already accreted onto the PNS. 
What is then seen in the equatorial wedge is the angular veloc- 
ity connected with material of initial cylindrical radii above 
~1000 km and, hence, of initially low angular velocity. How- 
ever, in the polar region we see material of large initial spher- 
ical radius, but relatively small (below ^1000 km) initial dis- 
tance from the rotation axis. Hence, the higher angular veloc- 
ity seen at 200 ms in the polar wedge. 

In the bottom-right panel of Fig. Owe display ml5b6:2D's 
angular velocity distribution at 200 ms after bounce. This 
model's initial central iron core period is quite long (31.73 s) 
and its iron core struct ure is similar to that of the older 
sll and sl5 progenitors (IWooslev & Weave3ll995h . Its iron 
core angular velocity profile exhibits the jumps at composi- 
tional boundaries which are characteristic of the set of ro- 
tating progenitors considered here. Again, arguing that the 
initial rotational profile is determining the PNS spin, we 
compare ml5b6:2D with si 1A1000P41.89:2D (bottom-right 
panel in Fig.llO>. Even though model ml5b6:2D has a shorter 
initial central period than si 1A1000P41.89:2D, its central 
spheroidal bulge of high-angular velocity is smaller at 200 ms 
after bounce. A glance at Fig. [^reveals the cause: ml5b6's 
iron core is indeed spinning faster than si 1A1000P41.89, 
though only in the inner ^700 km of the precollapse iron core. 
At that radius, ml5b6's O drops below the angular velocity in 
model si 1A1000P41.89:2D due to the use of the rotation law 
of eq. Q in the latter model. 

FigureFTSlshows equatorial angular velocity profiles at 200 
ms after bounce. The figure contains a panel with on a 
logarithmic scale out to 300 km in equatorial radius and a 
more zoomed-in version showing only the inner 100 km. ft 
is drawn on a linear scale to more clearly display the differ- 
ences between the individual models. We contrast the equa- 
torial profiles of models sllA1000:2D, sl5A1000:2D, and 
s20A1000:2D with Pq set to 2.34 s with the Heger models 
E15A:2D, E20A:2D, and ml5b6:2D. Above, we discussed 
the systematic changes of the PNS spin profile due to varia- 
tions in the initial rotational configuration for a fixed progeni- 
tor structure. In the discussion of the two-dimensional angular 
velocity snapshots, we saw that progenitor mass and iron core 
structure appear to have only a secondary effect on the PNS's 
birth spin rate and rotational configuration. This is supported 
by the equatorial rotational profiles in Fig. ^\ Models with 
similar initial (quantitative) rotational profile yield similar 



PNS rotational profiles, despite significant differences in their 
progenitor structures. Only in the linear-scale panels does one 
note differences that can be attributed to differences in the ini- 
tial angular velocity distribution ("spherical" vs. "cylindrical" 
and rotation law vs. one-dimensional rotation). For example, 
we noted the intere sting shape of E20A:2D's angular velocity 
distribution in Fig. 1141 The early drop of fl in the equato- 
rial regions is reflected in the dropoff in equatorial angular 
velocity at about 30 km (clearly visible in the right panel of 
Fig. I15> . Note that the initial discrete jumps of the angular 
velocity of models E15A:2D, E20A:2D, and ml5b6:2D are 
naturally smoothed out during the two-dimensional collapse 
and postbounce evolution. 

In Fig. [^_we compare the moment-of-inertia-weighted 
mean period Ppi2 of models with different progenitor struc- 
tures. For comparison, we enlarge the sample of models to in- 
clude sllA1000P10.47:2D and E15A:2Ds. The evolution of 
Ppi2 supports our assertion that progenitor structure has only 
a minor effect on the spin period of neutron stars at birth in 
our two-dimensional models. The magnitude of initial angu- 
lar velocity in the part of the precollapse iron core that forms 
the homologous inner core during collapse is the most im- 
portant parameter for the PNS spin. However, this assertion 
holds only for the early stages of the PNS evolution. It is pos- 
sible that the rotational configuration could be influenced by, 
among other phenomena, late-time fall-back accretion whose 
rate and absolute accreted mass will depend on the progenitor 
mass density profile. Furthermore, black hole formation will 
become more probable with increasing progenitor mass, lead- 
ing to a totally different scenario, possibly to a collapsar-type 
GRB (IWoosley|1993HMacFadyenet al.l2001h . 

To conclude this section on our 2D axisymmetric models, 
we discuss the evolution of the rotation parameter /3 for dif- 
ferent progenitor masses. In Fig. [n] we compare the post- 
bounce /3 evolution for a representative set of models. /3 is 
an integral quantity, evaluated by integrals that encompass 
all grid points. Differences (and similarities) in the /? evo- 
lution caused by progenitor structure and choice of initial an- 
gular velocity distribution and mapping are obvious. Models 
si 1 A1000P2.34:2D and sl5A1000P2.34:2D are evolved with 
the same initial rotation law and start with identical Pq. Their 
iron core structure is very similar as well. Their (3 behavior 
turns out to be almost identical. Model s20A1000P2.34:2D is 
evolved using the same initial rotational configuration, but has 
a significantly more massive iron core with an initially lower 
central density. Its higher value of (3 at bounce and during its 
postbounce evolution is a consequence of the much greater 
angular momentum contained in s20A1000P2.34:2D's core 
and of its relatively high initial value of f3. In Table |3] we 
list the total angular momentum contained in material above 
10'^ g cm~^ at the end of each model's numerical evolu- 
tion. The high value of (3 at bounce is also in part a con- 
sequence of its initially lower core density. This results in 
a greater spin up of the inner core material during collapse. 
Model E20A:2D's initial central spin period is a bit shorter 
than s20A1000P2.34:2D's (2 s vs. 2.34 s), while its initial in- 
ner core density is slightly larger E20A:2D's f3 is not much 
larger at bounce than s20A1000P2.34:2D's, but starts increas- 
ing with a greater rate postbounce. However, that increase is 
stopped about 10 ms after bounce when the initial profile's big 
jump in O at about 1000 km (Fig.[0 reaches small radii and 
the increase of rotational kinetic energy is abruptly slowed. 
Qualitatively, and to some extent quantitatively, this behav- 
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ior can be inferred from the two-dimensional angular velocity 
snapshot presented in Fig. 1 141 

Overall, E15A:2D is the fastest model of our two- 
dimensional model set. Its initial central period is as low as 
1.5 s and the initial value of /? amounts to 0.5%. Its initial 
inner core density is considerably lower than for the slower or 
nonrotating 15 Mq models (Fig.E). It does, Hke E20A;2D, 
exhibit a discontinuity in O in its initial profile. However, 
this discontinuity is located at greater radii and does not ad- 
vect in before the end of our simulation. All this explains the 
high /3 at bounce and the postbounce increase seen for this 
model. However, understanding what leads to the differences 
between models E15A:2D and E15A:2Ds requires a bit more 
involved discussion. In "shellular" rotation, the iron core an- 
gular velocity is set up to decrease roughly with r~^, where r 
is the spherical radius (see eq.[2and |3j- In comparison with 
"cylindrical" rotation, off-equatorial grid points at a given dis- 
tance !I! from the rotation axis are set up with an initial angular 
velocity that is reduced by a factor of (53/ \/05- + z-, where z is 
the cylindrical coordinate on the symmetry axis. Hence, there 
is less total angular momentum on E15A:2Ds' grid and its ini- 
tial value of /3 is about 30% smaller than E15A;2D's. These 
facts determine this model's /3 evolution. 

As expected, model E15A:2D reaches the greatest post- 
bounce /? of the entire set of two-dimensional models pre- 
sented here. At the end of the calculation, /3 is 12.85% and 
is still increasing. During their late-postbounce evolution, 
models E15A:2D and s20A1000P2.34:2D and, perhaps, mod- 
els E20A:2D and E15A:2Ds are likely to reach the classical 
threshold for secular rotational instability. However, it is un- 
likely that any of the models presented here will ever reach 
the threshold (3 at which classical high-T/| W| rotational insta- 
bilities would set in. Nevertheless, there are multiple models 
that may become unstable to the low-T/| W| instability of post- 
bounce cores observed by Ott et al. (2005). 

5. COMPARISON OF ONE- AND TWO-DIMENSIONAL 
SIMULATIONS 

In ^4. II and ^4.21 we presented the results of our one- and 
two-dimensional simulations separately. Here, we compare 
the two sets of simulations. Panel (a) of Fig. 1181 shows the 
equatorial angular velocity profiles at 200 ms after bounce 
for models si 1A1000P2.34, si 1A1000P10.47, and E15A. 
Two-dimensional variant (:2D) profiles are the dashed lines, 
while one-dimensional variants (:1D) are the solid lines. For 
model E15A we also plot the version with "shellular" rotation 
(:2Ds). First, we note that the rotational profiles of E15A:2D 
and E15A:2Ds lie practically on top of each other. This is 
remarkable, but not entirely surprising if one recalls that the 
initial equatorial angular velocity profiles of these models are 
identical. However, Fig. ^^demonstrates that there are signif- 
icant differences once one leaves the equator. 

Comparing the :1D profiles to their :2D counterparts, we 
note that there is good agreement outside ^50 km. With de- 
creasing radius the : ID and :2D profiles for a given model be- 
gin to diverge. The two-dimensional models have less smooth 
profiles. The substructure in the equatorial profiles is created 
by angular momentum advection in convective plumes inside 
the protoneutron star and in the postshock region. Since con- 
vection is confined to small cylindrical radii in quickly spin- 
ning PNSs, the equatorial profiles of faster models overall are 
smoother and closer to those of their :1D counterparts than 



those of slow models that allow for more vigorous convection 
at greater cylindrical radii. Note that the jump in fl visible in 
the ; ID profile of E15A (green solid-line, at ~20 km, also see 
Fig-GJ is smoothed out in its :2D counterpart. 

In panel (b) of Fig. we display the moment-of-inertia- 
weighted mean period Ppi2 for the same set of models. Ppi2 is 
an integral quantity based on a density criterion which takes 
into account the greater dimensionality of the :2D models. 
Model E15A is the fastest model for which we have :1D and 
:2D simulation data. Its initial central period is 1 .5 s and it un- 
dergoes core bounce under strong centrifugal influence. The 
postbounce Ppi2 evolutions of E15A:2D and E15A:2Ds dif- 
fer significantly from those of their one-dimensional counter- 
parts. At and after bounce, the two-dimensional PNSs are able 
to relax to greater PNS radii and energetically more favorable 
configurations than the spherically symmetric model. Hence, 
the :2D Ppi2 of this (almost) rotationally supported model is 
larger than the :lD's. With increasing initial rotation period, 
rotational support for the PNS becomes less and less and the 
overall greater angular momentum in the :2D models leads 
to shorter PNS Ppn s (see |3j "cylindrical" vs. " shell ular" 
rotation; also see the discussion towards the end of ^4.2t . 

The evolution of the rotation parameter f3 is shown in panel 
(c) of Fig. Again, there are significant differences be- 
tween :1D and :2D variants. For E15A, (3 is overestimated 
by the one-dimensional model and for the slower models 
S11A1000P2.34 and si 1A1000P10.47 the one-dimensional 
models underestimate [3. Model E15A:2Ds which is in "shel- 
lular" rotation reaches lower (3 than E15A:2D because there is 
less total angular momentum in the former The explanation 
for the differences in (3 between one-dimensional and two- 
dimensional models is analogous to the explanation for the 
differences seen in the P^n evolution discussed above. 

In panel (d) of Fig. ^] we display the mass- 
(spherical)radius relationship for selected :1D and :2D mod- 
els. Here, the naive first assertion holds: the slower a :2D 
model is, the more similar is its mass distribution to that of 
the corresponding :1D model and the better the :1D and :2D 
mass profiles agree. However, keep in mind that for fast :2D 
models that experience significant rotational flattening a map- 
ping M{r) will always be ambiguous. 

We conclude that the overall characteristics and systematics 
of the mapping between initial model and rotational profile to 
the postbounce spin of the PNS are rendered with acceptable 
accuracy by the one-dimensional treatment of rotation. There 
are, however, significant qualitative and quantitative differ- 
ences between :2D and :1D in all models. In particular, we 
do not expect convergence of the :1D and :2D results with 
increasing Po since multi-dimensional effects, especially con- 
vection that becomes more vigorous in slower models, will 
always distinguish :2D evolution from :1D. 

6. SIMPLE ESTIMATES 
OF THE FINAL NEUTRON STAR SPIN 

The simulations presented in this paper cover only the very 
early stages of a neutron star's life. During tens to hun- 
dreds of seconds after bounce and a successful explosion, 
the neutron star will deleptonize and cool via neutrino emis- 
sion. Fb£_the_typical neutron star mass range of ~ 1.2- 1. 6 
Mrr, (iFrver & Kalogera.200kiThorsett & Chaigabartvl 1999ft . 
iLattimer^R^kashl J2001l)' find a radius of ~12 km and a 
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fiducial moment of inertia of / = 0.35 MR^ for a cold neu- 
tron star Cold neutron stars are likely to be in solid-body 
rotation since this is the energetically most favorable configu- 
ration and viscous processes will most likely drive the neutron 
star into rigid rotation during its cooling phase. 

A number of studies of stellar evolution ("Heger et al.'200fll 
|2003b, 20 05) and of rotating core-collapse and supe rnova 
evolution dFrver & Hegeill2000t iFrver & Warrenll2004l) have 
estimated the final neutron star spin period via 



PNs = 2n- 



1 



(J /I) 



(5) 



Here J is the supposed final total angular momentum and / 
is the moment of inertia obtained with the above prescrip- 
tion and a choice of neutron star mass and radius (typically 
M = 1.4 M0 and R=12 km). For comparison with other 
studies, we estimate the asymptotic final neutron star spin of 
our models in the following manner. We consider our two- 
dimensional simulations and a representative subset of our 
one-dimensional simulations and choose M to be the mass of 
all material with a density above 10'^ g cm""* at 200 ms af- 
ter bounce. We choose this time to ensure that we have one- 
and two-dimensional data for our estimates. Also, that this is 
a sensible choice becomes apparent when looking at Fig. |8l 
which shows the mass interior to a given spherical radius at 
about 200 ms after bounce for representative models. Mpi2,2oo 
encompasses the majority of the total mass on the grid, lo- 
cated interior to some tens of kilometers. We assume that a 
successful explosion removes the remaining matter and that 
A^pi2.200 is the final neutron star mass. For simplicity, no fur- 
ther fall-back accretion or neutrino losses are considered. For 
the angular momentum in eq. (|5j, we choose Jpn.ioo, which is 
the total angular momentum of mass Mpi2.200- We assume this 
angular momentum and mass to be conserved during the con- 
traction to a final radius of 12 km. In addition, we compute 
the neutron star's rotation parameter /3ns = T/| W| and the ratio 
of the neutron star's angular velocity to the Keplerian break- 
up velocity, defined by ftK = y^GM/r^. For /3, we estimate 
the gravitational potential energy according to the prescrip- 
tion given in Lattimer & Prakash (2001 ): 



GM 



.GM, 



iWlc^Mc^x 0.6^/(1-0.5^) 



(6) 



We summarize our estimates in Table @] Neutron stars with 
/? > 50% (Viral theorem limit) and angular velocity above 
the Keplerian break-up velocity cannot exist. For six of 
the two-dimensional models listed in Table |4] our estimate 
leads to a spin rate beyond break-up and to unphysically 
high values of (3. Since the angular momentum contained in 
one-dimensional models is smaller than in the ;2D models, 
only the two most extreme :1D models reach such high (3 
and n I riK- These models simply cannot form stable compact 
neutron stars, unless there are one or multiple mechanisms 
that remove angular momentum from the cooling neutron 
star Otherwise, the neutron star will cool, but remain cen- 
trifugally hung-up with a disk of matter at Keplerian velocity. 
Two-dimensional models with initial central iron core spin 
period above ~5 s may form stable neutron stars below the 
Keplerian limit (si 1A1000P4.69:2D, si 1A1000P10.47:2D, 
si 1 AooP10.47:2D,sl 1 A1000P41.89:2D,sl 1A1000P157. 10:2D, 
and ml5b6:2D). For the one-dimensional models, this initial 
period cut is somewhat lower, allowing iron cores with 
initial periods of down to ^2 s to form stable neutron stars. 



However, keep in mind that the one-dimensional description 
of rotation is only approximate and underestimates the 
angular momentum in the PNS (see |5j- 

Only the overall slowest two-dimensional model 
(sllA1000P157.10;2D) produces a neutron star with a 
period in the ballpark of observed young pulsar periods 
(tens of milHseconds). Models si 1A1000P41.89:2D and 
ml5b6;2D yield neutron stars with ~7 ms period. This 
number is i n rough accord with the periods derived by 
iHeger et'all J2003bh for progenitors with similar initial 
rotation period and angular momentum distribution. How- 
ever, tHegeretal. (2005) estimated a neutron star period 
of 11 ms f or their ml5b6 prog enitor (cf. ml5b6:lD's P^s 
~ 14 ms). iFrver & WarrenI (12004) estimated for E15A and 
for a ml5b6-like model neutron star periods of 0.91 ms and 
17 ms, respectively. However, we note that in their three- 
dimensional simulations the PNS appears to continuously 
lose angular momentum after bounce. These authors do not 
provide a conclusive explanation for this phenomenon. 

We have surveyed the parameter space of initial central 
iron core periods from ~1.5 s to ~160 s with one- and two- 
dimensional supernova simulations. Across the board, our es- 
timates of the final neutron star spin periods, based on the 
two-dimensional simulations that handle rotation in a consis- 
tent manner, are shorter than found in previous studies. A 
large fraction of the two-dimensional models considered in 
this study produce cold, contracted neutron stars beyond the 
break-up spin rate. Given that observed young pulsar periods 
are between tens and hundreds of milliseconds, we conclude 
that there must be robust processes that efficiently spin down 
a neutron star before it emerges as a pulsar or that progenitor 
iron cores generically rotate more slowly than many of our 
progenitor models. We discuss in the following subsections a 
number of physical processes which could carry away or re- 
distribute a significant amount of the angular momentum of 
the nascent neutron star 



6.1. Viscous Processes 

In this paper, we have already considered a mechanism for 
spindown of PNSs: viscous angular momentum redistribution 
(by, e.g., the MRl)^. Because rotational energy scales with Q^, 
viscous heating scales with fl^, and the viscous timescale is 
proportional to il"', the multi ple effects of viscos ity are most 
noticeable in fast rotators (liThompson et alJl2005h . However, 
even for our rapidly rotating models (Pq 4 s), our results (see 
^4. H and Tables|2]and|4} indicate that viscous effects are not 
able to increase the spin period by more than a factor of 2 - 3 
on 100 ms to 1 s timescales. Although potentially interesting, 
these processes are not sufficient to yield fully-contracted 10- 
20 ms neutron stars. For slower initial iron core spin periods 
(Po Z 4 s), viscous effects quickly become much less impor- 
tant on the timescales considered here. 



6.2. Secular and Dynamical Rotational Instabilities 

PNSs that spin with high enoug h mig ht undergo the low- 
T/|W| instability (forT/|W| > 8%. IOtt et al. 2005) or the clas- 
sical secular and dynamical rotational instabilities at higher 

^ Other processes, such as the viscosity associated with convective tur- 
bulence, may also operate on the shear energy stored in differential rotation 
(Thompson et al. 2005) 
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T/|W|. These instabilities will put a lower limit on the neu- 
tron star spin of ~1 ms (set by the value of /3 at which the 
star regains stability, which in turn depends on the degree of 
differential rotation and the equation of state), but cannot lead 
to a significant spindown of the neutron star. 

6.3. r-Modes/Gravitational Wave Emission 

Rossby-waves (r-modes) in rapidly rotating neutron stars 
are driven unstable by the emission of gravitational radia- 
tion. Consequently, stellar rotational energ y is con verted into 
both gravitational wav es and r-mode energy jAnde rsson 1998; 
Lindblom et al. 2001). The spindown torque presented in 
Arras et al.i (.2003.) is 

7W ~ (2 X lO^^dyncm) 4],^ , (7) 

where is the saturation amplitude of the r-mode in dimen- 
sionless units and i'khz is the neutron star spin frequency (1/P) 
in units of kHz. With this, we estimate 



Bec ause of the low saturation level, '--^O.l JArras et all 
l2003i) . and the extreme spin-frequency dependence, r-modes 
can lead to a factor of two increase within a matter of days 
in the period of neutron stars that spin close to break-up, but 
will requi re millennia to slow down moderatel y rotating neu- 
tron stars ( lArras et al l2()()3tlHegeretalJ2005ft . 

6.4. Neutrino Emission 

Neutrinos carry away ~99% of the gravitational bind- 
ing energy of a forming neutron star Besides the possi- 
bility of a net linear momentum flux, which could be re- 
sponsible for the observed neutron star bir th kicks, neu- 
trinos also carry away angular momentum llEps tein"1978t 
iBaumgarte & Shapiro 19 98: Janka 2004. Assuming that the 
neutri nos are diffus ing out of the neutrinosphere in all direc- 
tions, Eankal (EqqI estimates that a maximum of ^40% of 
the total angular momentum can be lost if the axis ratio of the 
PNS is close to one. Hence, this mechanism could — in the 
most optimistic case — lead to a factor-of-two spindown for 
a relatively slow PNS (P^n > 10 ms). 

6.5. Rotation-Powered Siipernovae 

The typical supernova explosion kinetic energy is near 
10^' erg. The rotational energy stored in the PNSs in our mod- 
els is easily calculated using T = l/2yp.i2 (27r/Ppi2 ). Model 
s20A1000P2.34:2D has a rotational energy of ~2 x 10^- erg. 
If only 10% of this model's rotational energy is converted into 
radial kinetic energy, this could power the entire supernova 
explosi on, as suggested in the "jet-driven" supernova m ech- 
anism ( lAkivama et al ."^003': " Xkivama & Whee'ie3l2005l) . In 
that picture, the PNS represents the fly-wheel that provides 
the rotational energy to drive (magneto-) hydrodynamic jets 
along the poles leading to a bipolar supernova explosion. If 
10% of the rotational energy is dumped into the supernova, 
the PNS will spin down by not more than ~5%. If, however, 
50% of the rotational energy is imparted to the supernova, a 
spindown by ^^30% results. 



6.6. Magnetic Protoneutron Star Winds 

In the seconds after collapse and explosion, a thermal non- 
relativistic neutrino-driven wind emerges from the PNS (Dun- 
can et al. 1986). If the PNS has a magnetar-strength (> lO''*- 
10'^ G; Duncan & Thompson 1992) large-scale surface mag- 
netic field, then the wind material will be forced to coro- 
tate with the PNS. For rapid rotation, the wind is magneto- 
centrifugally whipped off of the stellar surface and the neutron 
star's angular momentum and rotational energy may be effi- 
ciently extracted (.Thomp_so_n e t al. 2004; Thompson 2003bj). 
The spin period increases as P/P w 3(M/M)(Ra/Rns)^, where 
M is the mass loss rate and Ra is the Alfven radius, where the 
kinetic energy density of the wind first exceeds the magnetic 
energy density. If the wind carries away angular momentum 
at a constant rate over a time t, then the relative change in pe- 
riod is P/Pbirth ~ exp(Pf/P). Assuming a 1.4 Mq PNS with a 
2 ms period, a radius of 20 km, an Alfven radius of 40 km, and 
a mass-loss rate of 10""' M0 s"' (jThomcson et al. 2004), to- 
gether with a 10 s wind duration, we obtain AP/P/„„/, of order 
10%. Although this mechanism does not spin the PNS down 
significantly on a 10 s timescale, the non-re lativistic wind car- 
ries of order 10^' erg in kinetic energy. 

As any PNS cools and contracts, the neutrino luminosity de- 
creases and the mass loss rate abates. As a consequence, for 
fixed magnetic field strength, the wind becomes Poynting-flux 
dominated and relativistic. For near-millisecond spin periods 
and magnetar-strength fields, this tra nsition happens just sec - 
onds after the explosion commences ([Thompson et alJ2004^ . 

Assuming a dipolar field configuration — probably pes- 
simistic considering the fact that pulsar braking indices are 
observed to be less than 3 (and, see Bucciantini et al. 2005) 
— the rotational energy loss rate is £ w B~R^il'^/c^ so that the 
spindown timescale is just P/P ~ 1.7 x lO'* s for a protoneu- 
tron star with R = 20 km, B=10''*G, P = 2 ms, and moment 
of inertia / = 0.35MR^. For these parameters, essentially all 
of the rotational energy can be lost to the relativistic wind in 
just tens of hours. 

The two most obvious and fundamental problems with spin- 
ning down a millisecond PNS rapidly with non-relativistic or 
relativistic pulsar-like winds is that (1) normal neutron stars 
are not expected to be born with large-scale magnetar-strength 
fields and (2) if the PNS is efficiently spun down, the remnant 
should bear the signature of the ^^10^" ergs of rotational en- 
ergy deposited within it. However, observations of generic 
young SNRs do not show evidence of such energy deposition. 

6.7. Fall-Back / Propeller Mechanism 

After a successful explosion there might be a phase dur- 
ing which the neutron star accretes a considerable amount of 
material that has failed to reach escape vel ocity in the shock- 
driven expulsion of the stellar envelope. iMacFadven et alJ 
(l2001h estimate the late-time (>1000 s after explosion; 
iHeger et'al]l2005l) fall-back accretion rate, based on the pa- 
rameters of SN 1987 A, to be 

M = 2xl0~^tf^\s-\ (9) 

where is the time in units of 10^ s. If the then deleptonized, 
but possibly centrifugally-hung-up neutron star has a suffi- 
ciently strong magnetic field, its Alfven radius can be greater 
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then the radius at which the corotation angular velocity ex- 
ceeds the Keplerian orbital speed. The infalling material is 
halted by magnetic forces at the Alfven radius, spun up to 
corotation and, if 17 > il^, expelled. The ejected matter car- 
ries away angular momentum and leads to a spindown of the 
neutro n star. iHeger et al. ( 20Q5). based on the work of lAlpad 
J2001h (also see Roman ova et al 1 120041 for recent numerical 
results), estimate that this so-called propeller mechanism can 
very efficiently spin down neutron stars with s urface mag- 
netic field strengths above ^2x10'^ G. Followin dHeger et al.l 
( l2005l) we estimate the integrated deceleration as 

AnK3Q.6fi%^nl^\ (10) 

where /i3o = B\2R\, B\2 is the surface field in units of 10'^ G, 
and /?6 is the neutron star radius in units of 10^ cm. O3 is the 
neutron star angular velocity in units of lO"' rad s"' . For a 2 ms 
period neutron star, hung-up with a 20 km radius and a conser- 
vative value of the surface magnetic field strength of 10'^ G, 
the integrated change in Vl is ^320 rad s"', corresponding to a 
change in the period of only 0.23 ms. With 15 times the mag- 
netic field strength and all other parameters kept fixed, the in- 
tegrated change in Q, would amount to ~2800 rad s~', leaving 
behind a neutron star with a period of ^18 ms. Given these 
numbers, there is great potential for the propeller mechanism 
to lead to a significant, order-of-magnitude spindown of the 
neutron star. However, it is not clear whether the propeller 
mechanism works in the way suggested in the literature. Fur- 
thermore, there is general ignorance concerning the magnetic 
field strength and topology in neutron stars. Also, there may 
not be much, or any, fallback for the majority of supernovae. 

7. SUMMARY AND DISCUSSION 

We have systematically investigated the mapping between 
presupernova (rotational) structure and the spin periods and 
rotational profiles of protoneutron stars. For this, we have per- 
formed a large set of rotating supernova core-collapse simula- 
tions in spherical- and axi-symmetry, using Newtonian grav- 
ity. For spherically-symmetric runs, we have employed the 
Boltzmann radiation-hydrodynamics code SESAME which 
has recently been upgraded to treat rotation in an approxi- 
mate way. In addition, viscous dissipation and angular mo- 
mentum redistribution have been included in a subset of our 
simulations. For models evolved in axisymmetry, we have 
made use of the multi-group flux-limited diffusion variant 
of the radiation-hydrodynamics code VULCAN/2D. For our 
simulations, we have taken 11, 15, and 20 Mq (sll, sl5, 
s20) presupernov a models from the stellar evolution study of 
iWooslev & Weaver (1995) and put them into rotation via a 
rotation law that assumes constant angular velocity either on 
cylindrical or spherical shells. In addition, we have performed 
simulations with t he progenitor models E15A and E 20A of 
IHeger et al.l ( l2000l) and ml5b6 of fleg eretal] ( l200l . These 
progenitor models include a one-dimensional prescription for 
rotational evolution. We have considered initial central iron 
core spin periods between 1.25 s and 157.1 s and only cen- 
tral periods that map to angular velocities below the Keplerian 
break-up limit at all radii in each model at all times during the 
numerical evolution. The initial range of the rotation param- 
eter /3 is set by a minimum of 7 x 10"^ % and a maximum of 
0.88%. 

In our one-dimensional simulations we find a roughly lin- 
ear relationship between initial central iron core spin period 



and the protoneutron star spin. As an analysis tool we de- 
fine a moment-of-inertia-weighted period (Ppi2 ) for all ma- 
terial above 10'^ g cm"^. In Fig.^] we plot Ppi2 as a func- 
tion of initial central spin period for a representative subset 
of our models. Each one-dimensional model (;1D) is marked 
with a solid box in the Pq - Ppi2 plane. The boxes of our 
sllA1000:lD model series are connected with solid lines to 
emphasize the Pq Ppi2 mapping. The graph is roughly fit- 
ted by a linear function with slope of 1.3 ms s~'. Only in the 
fastest models (initial periods shorter than ~2.5 s) do centrifu- 
gal effects lead to a noticable deviation from the displayed 
simple linear relationship. In considering the discrete profiles, 
we identify a region of near-solid body rotation inside the 
nuclear-density region of the PNS out to a radius of ^^10 km 
or ^0.5-0.7 M0 in mass coordinate. This is a consequence 
of quasi-homologous collapse. Outside this region PNSs are 
strongly differentially rotating. In a typical fast model, the 
angular velocity drops from ~2500 rad s~' at 20 km to a few 
hundreds rad s"' at 100 km. The energy stored in differential 
rotation can efficiently be tapped by viscous processes. We 
have considered the MRl as one such process and have in- 
cluded its action in a sub set o f our one-dimensional models 
(the ":lDv" models; see ^2.1> . As displayed in Fig.|5] vis- 
cosity leads to a spin down of the central part of the PNS and 
enforces rigid rotation out to 50 km and beyond. The viscous 
timescale is proportional to 17"' and the viscous heating rate 
scales with £7^. Hence, the efficiency of viscosity in reduc- 
ing the rotational shear, in moving angular momentum out, 
and in spinning down the PNS is only great for rapidly rotat- 
ing models. In Fig. open rhombi symbolize the models 
that include viscous angular momentum redistribution. Those 
rhombi corresponding to models of the sllA1000:lD series 
are connected with dashed lines. Comparing the two graphs 
in Fig. ^] we note that the inclusion of viscosity leads to 
slower Ppi2 for all models. The relative spin down is greatest 
for models below ~3 ms period, but no model is slowed by 
more than a factor of two in Ppi2 . Models with long initial 
periods are slowed by only a few to ten percent. 

We have compared the one-dimensional models with our 
two-dimensional, axisymmetric simulations that handle rota- 
tion consistently. We find that the ID approach correctly re- 
produces the overall phenomenology of the mapping between 
initial model and rotational profile to postbounce PNS spin. 
However, we identify non-negligible qualitative and quanti- 
tative differences when comparing individual models evolved 
in ID and 2D. In fast two-dimensional models convective mo- 
tions are confined to small cylindrical radii. Hence equatorial 
rotational profiles of fast 2D models exhibit less substructure 
and are closer to their ID counterparts than slower models in 
which convection is more vigorous out to larger radii. For 
slow to moderate initial spin rate, the two-dimensional sim- 
ulation results also yield a mapping of Pq Ppi2 which is 
roughly linear for large Pq and that deviates from linearity at 
lower values of Pq; the correct treatment of centrifugal effects 
results in less compact and more oblate PNSs. 

In comparing our 2D radiation-hydrodynamics, MGFLD 
models with purely hydrodynamic results, we find that 
the large-scale, coherent, and repeated expansion-collapse- 
bounce cycles observed in purely hydrodynamic calculations 
are critically damped in simulations that include neutrino ra- 
diative transfer and microphysics. This is an important finding 
because of its implication for the gravitational wave signa- 
ture of rotating stellar core collapse and bounce. Contrary to 
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curre nt belief ([Zw erger & Miilleiir i997t iDimmelmeier et alJ 
120021 iKotake et al.ll2003t lOtt et al.ll2004i) . the collapse and 
bounce of even fast rotators (like E15A with Po=1.5 s and 
an initial [3 of 0.42%) will not yield gravitational waves of 
large amplitude and energy emission. As with slowly rotating 
models, the greatest part of the gravitational wave emission 
will come from postbounce large-sc ale convective m otions in 
the postshock region and in the PNS dOtt et al.l2004h and, pos- 
sibly, from anisotropic n eutrino emission (iBurrows & HavesI 
[T996tlMuller et al.l2004 . 

At the end of our two-dimensional simulations (on average 
240 ms after core bounce), we find mean PNS periods in the 
range from 4.7 ms to 140 ms. The final P^n depends mainly 
on the inner precollapse iron core spin period. Differences in 
progenitor mass, structure, and rotation law lead to only first- 
order corrections in Ppi2 . However, differences in progenitor 
mass and structure have a strong impact on other observables; 
the more massive a progenitor, the greater are its final PNS 
mass, angular momentum, and rotation parameter (3. The 
maximum values of /3 range between 5xlO~^% and 12.9% 
in our two-dimensional simulations. At postbounce times, (3 
increases roughly linearly as the PNS simultaneously accretes 
more material and becomes more compact. During its cooling 
and deleptonization phase, a PNS with a birth period below 
ms will surpass the /3-threshold values for classical sec- 
ular and dynamical rotational instability (~14% and ~27%, 
respectively). It is, however, poss ible that they e xperience a 
low-T/|W| instability even earlier jOtt et alJ2005h . 

To compare with previous, though more approximate, stud- 
ies, we have taken the PNS configurations of our simulations 
at 200 ms after core bounce and extrapolated the spin periods 
and /3s to final, deleptonized, and cold states. For the same 
(or very similar) initial conditions, our m odels yield shor ter 
neutron star periods t han previous studies (Heaer et al. 200QI 
iFrver & Heger..2000t iHe ^er et al.ll2003b. .2005) . Any model 
with initial iron core period below '^lO s leads to a neu- 
tron star with a submillisecond period, spinning with surface 
velocities close to break-up. Since young pulsars enter the 
P - P dia gram at periods of man y tens to hundreds of mil- 
liseconds (iKaspi & Helfandl2002h . we have surveyed a num- 
ber of possible spindown mechanisms that could transform 
a young millisecond-period neutron star into one with a pe- 
riod in agreement with observations. These include viscous 
angular momentum redistribution by the MRl or another vis- 
cous process, rotational hydrodynamic instabilities, r-modes, 
neutrino angular momentum transport, rotation-powered ex- 
plosions, early magnetic-dipole (pulsar) emission, neutrino- 
driven magnetic winds, and the propeller mechanism. We find 
that none of the processes considered is robust enough to spin 
down a millisecond PNS to a tens-of-milliseconds pulsar The 
most promising is the propeller mechanism that removes an- 
gular momentum from the neutron star by expelling fall-back 
material at the Alfven radius, while not requiring magnetic 
field strengths in excess of several 10'^ G to a few 10'^ G. 
However, there remain many open questions concerning the 
physical details involved and whether fallback is even generic 
has not been resolved. 

Hence, we conclude that, if a pulsar were born rapidly ro- 
tating, nonaxisymmetric rotational hydrodynamic instabilities 
and other proprosed mechanisms of spindown may not be ad- 
equate to cause the requisite spindown from rapid to observed 
rates. In addition, it is clear that spindown mechanisms that 



employ magnetic fields would also transfer a rapidly rotating 
PNS's kinetic energy to the envelope debris, in contrast with 
the observational estima tes o f generic supernova and super- 
nova remnant energies ( ^6.6> . Therefore, on the basis of our 
theoretical calculations of the mapping between initial and fi- 
nal spins, likely spindown mechanisms, and the observational 
constraints from pulsars and supernova energetics, we con- 
clude that the progenitor spin rates of the cores of the massive 
stars that give birth to observed neutron stars and pulsars are 
low, with corresponding periods on average greater than ^50 
seconds. The associated average neutron star would thereby 
be born with a period no less than ~10 ms and with /3 < 2%. 
However, the final neutron star spin is a function of the pro- 
genitor mass. Some very massive progenitors might lead to 
very fast and massive neutron stars that might evolve to col- 
lapsars and rotation-powered GRBs. 

Figure |20| shows equatorial specific angular momentum 
profiles at 200 ms after core bounce for a representative sub- 
set of our models evolved in axisymmetry with VULCAN/2D. 
Note that models of similar initial central iron core Pq yield 
very similar specific angular velocity profiles, despite dif- 
ferences in progenitor structure. Model si lAooP10.47:2D, 
which is initially set up in solid-body rotation, yields out- 
side some tens of kilometers specific angular momenta as 
high as found in differentially rotating models with up to 
5 times shorter initia l cent ral periods. Given the estimates 
of Heaer & Woosley (l2003h . the specific angular momenta of 
our models with initial central periods below ~10 s are suf- 
ficient to rotationally support an accretion disk around a Kerr 
black hole. 

Finally, we point out that even though this study has 
advanced our knowledge concerning the birth spin peri- 
ods and rotational profiles of neutron stars, and their de- 
pendence on initial iron core spin and structure, many 
important questions remain unanswered. We have fol- 
lowed the PNSs for only a few hundred milliseconds af- 
ter formation. We have not tracked the supernova explo- 
sion through the stellar envelope and subsequent fall-back 
accretion. We have not included general relativistic ef- 
fects, nor magneto-hydrodynamics. Fully consistent, three- 
dimensional, radiation-magneto-hydrodynamics simulations 
in general relativity will be required to reveal in detail the 
evolution from a very young PNS to the final rotating pulsar 
(or black hole). 



We acknowledge discussions with and help from Jeremiah 
Murphy, Itamar Lichtenstadt, Casey Meakin, and Ed Sei- 
del. Importantly, we acknowledge support for this work 
from the Scientific Discovery through Advanced Comput- 
ing (SciDAC) program of the DOE, grant number DE-FC02- 
01ER41 184, and from the NSF grant number AST-0504947. 
E. L. acknowledges support from the Israel Science Founda- 
tion under grant 805/04. R.W. thanks the Institute of Astron- 
omy of the ETH Zurich for providing part-time office space. 
T.A.T. acknowledges support by NASA through Hubble Fel- 
lowship grant #HST-HF-01 157.01 -A, awarded by the Space 
Telescope Science Institute, which is operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., for 
NASA, under contract NAS 5-26555. Finally, we thank Jeff 
Fookson and Neal Lauver of the Steward Computer Support 
Group for their invaluable help with the local Beowulf clus- 
ters and acknowledge the use of the NERSC/LBNL/seaborg 



17 



and ORNL/CCS/RAM machines. Movies and additional still from the authors, 
frames associated with this work can be obtained upon request 



REFERENCES 



Akiyama, S., Wheeler, J. C, Meier, D. L., & Lichtenstadt, I. 2003, ApJ, 584, 
954 

Akiyama, S., & Wlieeler, J. C. 2005, ApJ, 629, 414 
Alpar, M. A. 2001, ApJ, 554, 1245 
Andersson, N. 1998, ApJ, 502, 708 

Arras, R, Flanagan, E. E., Morsink, S. M., Schenk, A. K., Teukolsky, S. A., 

& Wasserman, I. 2003, ApJ, 591, 1 129 
Balbus, S. A., & Hawley, J. R 1991, ApJ, 376, 214 
Balbus, S. A., & Hawley, J. R 1994, MNRAS, 266, 769 
Baumgarte, T. W., & Shapiro, S. L. 1998, ApJ, 504, 431 
Bethe, H. A., & Wilson, J. R. 1985, ApJ, 295, 14 
Bruenn, S. W. 1985, ApJS, 58, 771 

Bucciantini, N., Thompson, T. A., Arons, J., Quataert, E., & Del Zanna, L. 

2005, submitted to MNRAS 
Buras, R., Rampp, M., Janka, H.-T., & Kifonidis, K. 2003, Phys. Rev. Lett., 

90, 241101 

Buras, R., Rampp, M., Janka, H. -., & Kifonidis, K. 2005, ApJ. accepted, 

astro-ph/0507135 
Burrows, A., & Hayes, J. 1996, Phys. Rev. Lett., 76, 352 
Burrows, A., Hayes, J., & Rryxell, B. A. 1995, ApJ, 450, 830 
Burrows, A., & Lattimer, J. M. 1986, ApJ, 307, 178 

Burrows, A., Young, T., Pinto, P., Eastman, R., & Thompson, T. A. 2000, 
ApJ, 539, 865 

Centrella, J., New, K., Lowe, L., & Brown, J. 2001, ApJ, L550 

Dimmelmeier, H., Ront, J., & Muller, E. 2002, A&A, 393, 523 

Duncan, R. C, Shapiro, S. L., & Wasserman, L 1986, ApJ, 309, 141 

Duncan, R. C, & Thompson, C. 1992, ApJ, 392, L9 

Epstein, R. 1978, ApJ, 223, 1037 

Fryer, C. L. 2004, ApJ, 601, LI 75 

Fryer, C. L. 1999, ApJ, 522, 413 

Fryer, C. L., & Heger, A. 2000, ApJ, 541, 1033 

Fryer, C. L., & Kalogera, V. 2001, ApJ, 554, 548 

Fryer, C. L., & Warren, M. S. 2004, ApJ, 601, 391 

Fukuda, I. 1982, PASP, 94, 271 

Heger, A., Langer, N., & Woosley, S. E. 2000, ApJ, 528, 368 

Heger, A., & Woosley, S. E. 2003, AIP Conf. Proc. 662: Gamma-Ray Burst 

and Afterglow Astronomy 2001 : A Workshop Celebrating the First Year of 

the HETE Mission, 662, 214 
Heger, A., Fryer, C. L., Woosley, S. E., Langer, N., & Hartmann, D. H. 2003a, 

ApJ, 591, 288 

Heger, A., Woosley, S. E., Langer, N., & Spruit, H. C. 2003b, Stellar 

Rotation. Proceedings lAU Symposium No. 215 
Heger, A., Woosley, S. E., & Spruit, H. C. 2005, ApJ, 626, 350 
Herant, M., Benz, W., Hix, W. R., Riyer, C. L., & Colgate, S. A. 1994, ApJ, 

435, 339 

Ho, W. C. G., & Lai, D. 2000, ApJ, 543, 386 

Janka, H.-T., Buras, R., & Rampp, M. 2003, Nuc. Phys. A, 718, 269 
Janka, H.-T. 2004, in "Young Neutron Stars and Their Environments", 

Proceedings, lAU Symposium 218 
Janka, H.-T, Scheck, L., Kifonidis, K., Muller, E., & Plewa, T. 2005, in: The 

Fate of the Most Massive Stai's, Proc. Eta Carinae Science Symposium 

(Jackson Hole, May 2004) 
Kaspi, V. M., & Helfand, D. J. 2002, Neutron Stars in Supernova Remnants, 

ASP Conference Series, Vol. 271, edited by P O. Slane & B. M. Gaensler 
Keil, W., & Janka, H.-T. 1995, A&A, 296, 145 
Kotake, K., Yamada, S., & Sato, K. 2003, Phys. Rev. D, 68, 044023 



Lai, D., Chernoff, D. R, & Cordes, J. M. 2001, ApJ, 549, 1111 
Lattimer, J. M., & Prakash, M. 2001, ApJ, 550, 426 
Lattimer, J. M., & Swesty, R D. 1991, Nucl. Phys. A, 535, 331 
Liebendorfer, M., Mezzacappa, A., & Thielemann, R. 2001, Phys. Rev. D, 63, 
104003 

Lindblom, L., Tohline, J. E., & Valhsneri, M. 2001, Phys. Rev. Lett., 86, 1 152 
Livne, E. 1993, ApJ, 412, 634 

Livne, E., Burrows, A., Walder, R., Lichtenstadt, I., & Thompson, T. A. 2004, 
ApJ, 609, 277 

Lyne, A. G., & Graham-Smith, R. 1998, Pulsar astronomy, Cambridge, U.K.: 

Cambiidge University Press. 
Lyne, A. G., & Lorimer, D. R. 1994, Nature, 369, 127 
MacRadyen, A. I., Woosley, S. R., & Heger, A. 2001, ApJ, 550, 410 
Muller, E., Rampp, M., Buras, R., Janka, H.-T, & Shoemaker, D. H. 2004, 

ApJ, 603, 221 

Mai-shall, R E., Gotthelf, E. V., Zhang, W., Middleditch, J., & Wang, Q. D. 

1998, ApJ, 499, LI 79 
Martinez-Pinedo, G., Liebendoerfer, M., & Prekers, D. 2004, 
jastro-ph/04 12091 

Mezzacappa, A., Liebendorfer, M., Messer, O. E., Hix, W. R., Thielemann, 

R, & Bruenn, S. W. 2001, Phys. Rev. Lett., 86, 1935 
Monchmeyer, R., Schafer, G., Muller, E., & Kates, R. 1991, A&A, 246, 417 
Nomoto, K., Tominaga, N., Umeda, H., Maeda, K., Ohkubo, T, Deng, J., & 

Mazzali, P A. 2005, ASP Conf Ser. 332: The Rate of the Most Massive 

Stars, 332, 384 
Ostriker, J. P, & Gunn, J. E. 1969, ApJ, 157, 1395 
Ott, C. D., Burrows, A., Livne, E., & Walder, R. 2004, ApJ, 600, 834 
Ott, C. D., Ou, S., Tohline, J. E., & Burrows, A. 2005, ApJ, 625, LI 19 
Pons, J. A., Reddy, S., Prakash, M., Lattimer, J. M., & Miralles, J. A. 1999, 

ApJ, 513, 780 
Rampp, M., & Janka, H.-T. 2002, A&A, 396, 361 

Romanova, M. M., Ustyugova, G. V., Koldoba, A. V., & Lovelace, R. V. E. 

2004, ApJ, 616, L151 

Saijo, M., Baumgarte, T. W., & Shapiro, S. L. 2003, ApJ, 595, 352 

Scheck, L., Plewa, T, Janka, H.-T, Kifonidis, K., & Muller, E. 2004, Phys. 

Rev. Lett., 92, 011103 
Shibata, M., Karino, S., & Eriguchi, Y. 2002, MNRAS, 334, L27 
— . 2003, MNRAS, 343, 619 

Spruit, H. C, & Phinney, R. S. 1998, Nature, 393, 139 

Tassoul, J.-L. 1978, Theory of Rotating Stars (Princeton U. S. A.: Princeton 

University Press) 
Thompson, T. A. 2003a, ApJ, 585, L33 
— . 2003b, ApJ, 585, L33 

Thompson, T. A., Burrows, A., & Pinto, P A. 2003, ApJ, 592, 434 
Thompson, T. A., Chang, P, & Quataert, E. 2004, ApJ, 611, 380 
Thompson, T. A., Quataert, R., & Burrows, A. 2005, ApJ, 620, 861 
Thorsett, S. E., & Chakrabarty, D. 1999, ApJ, 512, 288 
Walder, R., Burrows, A., Ott, C. D., Livne, E., Lichtenstadt, I., & Jan-ah, M. 

2005, ApL 626,317 

Watts, A. L., & Andersson, N. 2002, MNRAS, 333, 943 
Woosley, S. E. 1993, ApJ, 405, 273 
Woosley, S. R., & Weaver, T. A. 1995, ApJS, 101, 181 
Zwerger, T, & Muller, R. 1997, A&A, 320, 209 



18 



TABLE 1 
Initial Models 



Progenitor 


ZAMS Mass 




Mfe 


Model Name 


Evolution 


A 





Pa 




(Mq) 


(10** cm) 


(Mq) 




Type 


(km) 


(rad s ) 


(s) 


sll 




1.35 


1 32 




















ail i \J\f\jL 1 .Z. J . 1 U 


Oijkj^^iVl±j 


1000 


5 03 


1.25 










« 1 1 A 1 nonp? no- 1 






J. it 


9 on 














1 ono 




9 14 










si 1 A1 000P3 on- in" 


SFSAMF 


1000 


2.09 


3.00 










sllA1000P4.00:lD'' 


SESAME 


1000 


1.57 


4.00 










sllA1000P5.00:lD" 


SESAME 


1000 


1.26 


5.00 










sllAlOOOPS.OOilD" 


SESAME 


1000 


0.79 


8.00 










sllA1000P10.47:lD 


SESAME 


1000 


0.60 


10.47 










sllA1000P2.00:lDv" 


SESAME V 


1000 


3.14 


2.00 










sllA1000P2.34:lDv 


SESAME V 


1000 


2.68 


2.34 










sllA1000P3.00:lDv" 


SESAME V 


1000 


2.09 


3.00 










sllA1000P4.00:lDv" 


SESAME V 


1000 


1.57 


4.00 










sllA1000P10.47:lDv 


SESAME V 


1000 


0.60 


10.47 










sllAooP10.47:lD 


SESAME 


oo 


0.60 


10.47 










sllAooP10.47:lDv 


SESAME V 


oo 


0.60 


10.47 










sllA1000P2.34:2D'' 


VULCAN 


1000 


2.68 


2.34 










sllA1000P2.34:2Dh 


VULCAN H 


1000 


2.68 


2.34 










sllA1000P4.69:2D'' 


VULCAN 


1000 


1.34 


4.69 










s11A1000P10.47:2D'' 


VULCAN 


1000 


0.60 


10.47 










sllA1000P10.47:2Dh 


VULCAN H 


1000 


0.60 


10.47 










s11A1000P41.89:2D'' 


VULCAN 


1000 


0.15 


41.89 










t;! 1 A 1 OOOPl S7 UV9r)t' 
ail 1 \j\jijr u / . i\j . 


VTTT CAN 


1000 


0.04 


157 10 










q1 1 A^PI n 47-9n 


VTTT PAN 


oo 


60 


10 47 


si S 


1 s 


I.ID 


1 28 




















sl5A1000P2.34:lD 


SESAME 


1000 


2.68 


2.34 










sl5A1000P10.47:lD 


SESAME 


1000 


0.60 


10.47 










sl5A1000P2.34:lDv 


SESAME V 


1000 


2.68 


2.34 










sl5A1000P10.47:lDv 


SESAME V 


1000 


0.60 


10.47 












VTTT PAN 


1000 


2 68 


9 H 


s20 


20 


2.21 


1.92 




















s20A1000P2.34:lD 


SESAME 


1000 


2.68 


2.34 










s20A1000P10.47:lD 


SESAME 


1000 


0.60 


10.47 










«9nA 1 nOOP? ^4- 1 Dv 

aZ.\Jrt. i \J\J\JL L }'-r. 1 \J V 




1000 


2 68 


2.34 










s20A1000P10.47:lDv 


SESAME V 


1000 


0.60 


10.47 










s20A1000P2.34:2D 


VULCAN 


1000 


2.68 


2.34 


E15A 


15 


2.22 


1.62 


E15A:1D'' 
E15A:lDv" 
E15A:2D 
E15A:2Ds 


SESAME 
SESAME V 
VULCAN 
VULCAN S 




4.18 
4.18 
4.18 
4.18 


1.50 
1.50 
1.50 
1.50 


E20A 


20 


2.78 


1.84 


E20A:1D 

E20A:lDv 

E20A:2D 


SESAME 
SESAME V 
VULCAN 




3.13 
3.13 
3.13 


2.00 
2.00 
2.00 


ml5b6 


15 


1.95 


1.47 


ml5b6:lD 

ml5b6:lDv 

ml5b6:2D 


SESAME 
SESAME V 
VULCAN 




0.198 
0.198 
0.198 


31.73 
31.73 
31.73 



Note. — List of models included in this st udy. The sup ernova progenitor star models sll, sl5 and s20 are taken from lWooslev & WeaveJ 119951) ; E15A and E20A 
are from lHeger erani200Cll and mI5b6 is from lHeger eTaiTlaOOl . Models E15A, E20A, and ml5b6 include rotation using a one-dimensional prescription for angular 
momentum evolution. Centrifugal force terms were included up to the end of core carbon burning, sll, si 5, and s20 are non-rotating and are forced to rotate by the rotation 
law given in eq. 0. Rfe is the radius of the iron core (determined by the discontinuity in Y^, at the outer edge of the iron core). Mfe is the mass of the iron core. Evolution type 
denotes the way a model is evolved. SESAME stands for spherically symmetric evolution with the SESAME code, SESAME V indicates that dissipative angular momentum 



redistribution is switched on, VULCAN means two-dimensional simulation with MGFLD, VULCAN S stands for two-dimensional simulation with MGFLD, but shelliilar 
rotation ( jHl . and VULCAN H means purely hydrodynamic two-dimensional simulation without neutrinos. A is the parameter governing the degree of differential rotation 
(eq.0. and Pq are the initial central angular velocity and spin period, respectively. All progenitor models already have an initial infall velocity profile when they are 
mapped onto our computational grid. Nothing aitificial is done to initiate collapse. 

^Previously published/discussed in lThompson et all 120051) 
''Previously published/discussed in lWalder et al J |20Q^ 
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TABLE 2 



Summary of Results : One-Dimensional Models 


Model Name 


t/ - tfc 


Po 






Jnll.f 


A- 


ft- 




(ms) 


(s) 


(ms) 


(Mq) 


(Kr' erg s) 


(%) 


(%) 




587 


1.25 


3.11 


1.07 


1.172 


8.80-10"' 


18.44 


s 1 1 A 1 onop? on- 1 n 


629 


2.00 


3.50 


1.28 


1.127 


3.44. 10"' 


10.45 


sllA1000P2.34:lD 


630 


2.34 


3.85 


1.32 


1.091 


2.51-10-' 


7.96 


sllA1000P3.00:lD 


588 


3.00 


4.72 


1.35 


0.901 


1.53-10-' 


5.10 


sllA1000P4.00:lD 


541 


4.00 


5.94 


1.37 


0.706 


8.60-10-2 


2.98 


sllA1000P5.00:lD 


608 


5.00 


6.94 


1.39 


0.586 


5.50-10-2 


1.97 


sllA1000P8.00:lD 


532 


8.00 


11.16 


1.39 


0.366 


2.15-10-2 


0.759 


sllA1000P10.47:lD 


549 


10.47 


14.28 


1.39 


0.283 


1.25-10-2 


0.448 


sllA1000P2.00:lDv 


500 


2.00 


6.72 


1.13 


0.392 


3.44- 10-' 


2.25 


si lA100()P2.34:lDv 


437 


2.34 


7.18 


1.16 


0.428 


2.5 l-K)-' 


2.42 


Si lAIUUUri.UU. lUV 


1037 


3.00 


6.13 


1.29 


0.500 


1 .53- K)-' 


1.97 


!■ 1 1 A 1 nnAD/l An- 1 t^u 
SI IAIUUUr4. UU. lUV 


J / J 


A nn 
4.UU 


/.42 


1 n 
I .J J 


U. J24 


c An 1 n— 2 
o.oU- lU 


1 OQ 


f 1 1 A 1 nnAD in An • \ t\,7 
Si lAIUUUr IU.4 /. IL)V 


549 


10.47 


14.93 


1.39 


0.271 


1 .25- 10-2 


0.416 


sllAooP10.47:lD 


571 


10.47 


9.18 


1.34 


0.428 


5.51-10-2 


1.59 


sllAooP10.47:lDv 


570 


10.47 


9.41 


1.34 


0.415 


5.51-10-2 


1.17 


sl5A1000P2.34:lD 


612 


2.34 


3.95 


1.39 


1.197 


2.45- 10-' 


7.92 


sl5A1000P10.47:lD 


562 


10.47 


13.70 


1.47 


0.316 


1.22-10-2 


0.468 


s 1 5 A 1 000P2. 34 : 1 D v 


612 


2.34 


6.76 


1.24 


0.501 


2.45- 10-' 


2.51 






1U.4 / 


1/1 11 

14. 1 J 


1 AA 
1 .40 




1 1 n^2 


n AAi 
U.44 / 


s20A1000P2.34:lD 


462 


2.34 


3.48 


1.84 


2.820 


3.53- 10-' 


15.26 


s2()A100()P10.47:lD 


349 


10.47 


11.30 


1.85 


0.645 


1.76-10-2 


0.875 


s20A1000P2.34:lDv 


365 


2.34 


6.80 


1.40 


0.749 


3.53- K)-' 


4.32 


s20A1000P10.47:lDv 


347 


10.47 


11.86 


1.85 


0.620 


1.76-10-2 


0.740 


E15A:1D 


780 


1.50 


2.91 


1.53 


2.064 


4.95- 10-' 


17.45 


E15A:lDv 


639 


1.50 


5.76 


1.25 


0.513 


4.95- 10-' 


2.50 


E20A:1D 


461 


2.00 


3.43 


1.84 


2.213 


2.85- K)-' 


12.52 


E20A:lDv 


361 


2.00 


6.92 


1.46 


0.805 


2.85- 10-' 


3.95 


ml5b6:lD 


589 


31.73 


57.98 


1.50 


0.077 


7.77-10^ 


0.0280 


ml5b6:lDv 


589 


31.73 


58.21 


1.50 


0.076 


7.77- 10-* 


0.0234 


Note. — Numerical results for the supernova simulations performed with SESAME in 


spherical symmetry. Model initial conditions 


are detailed in Table [Tl Here, tf - t/^ 


denotes the time elapsed after 


core bounce to when 


we stop the calculation. Po is the initial central rotation period and Pp\2.f is the mean PNS rotation period (for 


material with 


p > lO'" g cm ^ (defined by eq.^ at the end of the simulation. 


Mpi2,f is the integrated 


mass of all material above the same density threshold at the end of the numerical 


calculation and Jp\2.f' is the angular momentum associated with Mp 


2./. /3i and represent the values of T/|W at the initial stage and at the end of the simulation, 


respectively. 








TABLE 3 












Summary of Results : Two-Dimensional Models 






Model Name 


tf-tb 


Po 


Ppl2./ 


Mpn.f 


^12./ 


ft 


ft' 




(ms) 


(s) 


(ms) 


(Me) 


(lO"*' ergs) 


(%) 


(%) 


sllA1000P2.34:2D 


197 


2.34 


4.74 


1.17 


1.522 


3.14-10-' 


6.72 


sllA1000P2.34:2Dh 


150 


2.34 


4.86 


1.16 


1.436 


3. 14-10-' 


6.12 


sllA1000P4.69:2D 


262 


4.69 


5.63 


1.28 


0.998 


8.35-10-2 


2.64 


sllA1000P10.47:2D 


293 


10.47 


10.05 


1.33 


0.489 


1.53-10-2 


0.694 


sllA1000P10.47:2Dh 


294 


10.47 


10.44 


1.13 


0.301 


1.53-10-2 


0.498 


sllA1000P41.89:2D 


319 


41.89 


37.79 


1.33 


0.123 


1.04-10-3 


0.0536 


sllA1000P157.10:2D 


318 


157.10 


139.63 


1.33 


0.033 


7.00-10-5 


0.00451 


sllAooP10.47:2D 


192 


10.47 


9.50 


1.26 


0.613 


5.35-10-2 


1.26 


sl5A1000P2.34:2D 


215 


2.34 


4.69 


1.18 


1.563 


3.07- K)-' 


7.00 


s20A1000P2.34:2D 


238 


2.34 


4.98 


1.42 


3.351 


5.08- K)-' 


12.67 


E15A:2D 


217 


1.50 


5.12 


1.28 


2.801 


6.45- K)-' 


12.85 


E15A:2Ds 


243 


1.50 


4.73 


1.37 


2.738 


4.64- K)-' 


10.75 


E20A:2D 


211 


2.00 


5.40 


1.49 


3.185 


4. 17-10-' 


10.76 


ml5b6:2D 


195 


31.73 


50.25 


1.39 


0.129 


1.03-10-3 


0.0533 



Note. — 



Same as Tablel2l but listing results of our simulations performed with VULCAN/2D in axisymmetry. 
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TABLE 4 

FrNAL Neutron Star Spin Estimates 



Model Name 


Po 




Jpn.200 




13m 






(s) 


(M©) 


(10*' erg s) 


(ms) 


(%) 






2 34 


1 17 


1 522 


48 


50 65 


1 37 


sllA1000P4.69:2D 


4.69 


1.25 


0.992 


0.85 


15.16 


0.75 


si 1 A lonopi n 47-2n 


10.47 


1.28 


0.443 


1.82 


3.25 


0.35 


si 1A1000P41.89:2D 


41.89 


1.29 


0.113 


7.19 


0.21 


0.09 


si 1A1000P157.10:2D 


157.10 


1.29 


0.030 


27.10 


0.02 


0.02 


sllAooP10.47:2D 


10.47 


1.26 


0.613 


1.30 


6.54 


0.49 


sl5A1000P2.34:2D 


2.34 


1.18 


1.550 


0.48 


51.25 


1.38 


s2nAinnoP7 34-2n 


2.34 


1.27 


3.098 


0.28 


129.01 


2.20 




1.50 


1.27 


2.724 


0.29 


126.05 


2.17 


Fl^A-7r)s 


1.50 


1.33 


2.553 


0.33 


96.01 


1.89 


E20A:2D 


2.00 


1.47 


3.104 


0.30 


104.13 


1.98 


ml5b6:2D 


31.73 


1.39 


0.129 


6.79 


0.21 


0.09 


sllA1000P2.34:lD 


2.34 


1.15 


0.762 


0.95 


13.39 


0.70 


sllA100()P2.34:lDv 


2.34 


1.06 


0.378 


1.77 


4.23 


0.39 


sllA1000P10.47:lD 


10.47 


1.26 


0.211 


3.76 


0.77 


0.17 


sllA1000P10.47:lDv 


10.47 


1.25 


0.200 


3.94 


0.71 


0.16 


sllAooP10.47:lD 


10.47 


1.23 


0.295 


2.63 


1.63 


0.25 


sllAooP10.47:lDv 


10.47 


1.22 


0.275 


2.80 


1.45 


0.23 


sl5A1000P2.34:lD 


2.34 


1.19 


0.790 


0.95 


12.96 


0.69 


sl5A1000P2.34:lDv 


2.34 


1.09 


0.416 


1.65 


4.71 


0.42 


s20A1000P2.34:lD 


2.34 


1.41 


1.622 


0.55 


32.35 


1.10 


s20A1000P2.34:lDv 


2.34 


1.15 


0.527 


1.38 


6.41 


0.49 


E15A:1D 


1.5 


1.21 


1.271 


0.60 


31.86 


1.09 


E15A:lDv 


1.5 


1.01 


0.378 


1.68 


4.91 


0.42 


E20A:1D 


2.00 


1.49 


1.465 


0.64 


22.24 


0.92 


E20A:lDv 


2.00 


1.19 


0.560 


1.34 


6.51 


0.49 


ml5b6:lD 


31.37 


1.33 


0.059 


14.21 


0.05 


0.04 


ml5b6:lDv 


31.37 


1.33 


0.058 


14.45 


0.05 


0.04 



Note. — Overview table listing the key quantities for the cold neutron star spin estimates of ^ Pq 
is the initial central rotation period, Mp]2.200 is the mass of all material with density p > lO'" g cm~"^ at 
200 ms after bounce, and 7p 12.200 is ^he total angular momentum contained in that material. P^s is the final 
neutron star spin estimate according to eq. jsj. jS^^js is an estimate of the neutron star rotation parameter 
and ilj^s/^K is the ratio of the neutron star's angular velocity to the Keplerian break-up velocity defined 
by ^J GM/r^. Neutron stars with 0,^5 /^k > 1.0 or /9 > 50% cannot exist. 
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Initial Q Profiles 

E15A 

E20A 

ml5b6 

rotation law, A^IOOO kir 




\ b 



Fig. 1 . — Initial angular velocity of the rotating progenitor models E15A (magenta), E20A (orange), and ml5b6 (black) and of models that are set into rotation 
according to the rotation law of eq. Q, with A fixed to 1000 km (shown in green for various Po or f2o). The discontinuities in the E15A, E20A, and ml5b6 
profiles are comiected to compositional interfaces in the rotating progenitors (see .Heger et aL200Q,) . 




r (cm) 



Fig. 2. — Initial mass density profiles for all progenitor models used in this study, sll, sl5, s20 (II M p,, 15 Mq, and 20 Mq, r espectively) are taken from 
IWooslev & Weaver! H99l . Models E15A and E20A (15 M0and 20 M©) are taken fron JHeger eTai]i200Cl) and model ml5b6 is from lHeger eTalll200l . E15A, 
E20A, and ml5b6 were evolved in spherical symmetry with an approximate treatment of centrifugal effects until the end of core carbon burning. Note the great 
similarity in the density stratifications of models sll, sl5, and ml5b6, especially interior to ~2000 km. 
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I — ' — ' — ' — ' — I — ' — ' — ' — ' — I — ' — ' — ' — ' — I — ' — ' — ' — ' — I — ' — ' — ' — ' — ' 
100 200 300 400 500 

t - (ms) 

Fig. 3. — Evolution of the moment of inertia-weighted mean period Fpi2 for material of density greater than or equal to lO'^ g cm"' according to eqs. l3l 
and |2)- Shown are selected models of the one-dimensional sll model series with (dotted Hnes, :lDv) and without (solid lines, :1D) viscous dissipation/angular 
mom entum redistribution. For models si 1A1000P1.25:1D and si 1A1000P8.00:1D, no computation with active viscous dissipation is performed. Viscosity (see 
i|2.1l genetically leads to spindown of the protoneutron star. It is most efficient at doing so in fast rotators and in the most extreme case (in model sllA1000P2.00) 
leads to a final mean period that is twice as large as in the nonviscous case. 




— 9UA1000P1.25 — S11A1000P8.00 

— sUAlOOOPa.OO — sUAlOOOPlO.47 

— si 1A1000P3 00 — sllA«'P10.47 

— S11A1000P4.00 




200 400 

t - (ma) 



Fig. 4. — Evolution of the rotation parameter /3 = T/|W| for selected models of the sll model series in ID. Solid lines: models without viscous dissipation (:1D), 
dotted lines: models with viscous dissipation (:lDv). For models sllA1000P1.25:lD and sllA1000P8.00:lD, no computation with active viscous dissipation is 
performed. Note the strong impact of viscous dissipation acting on the energy stored in differential rotation. 
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Fig. 5. — Profiles of the period (a) and the angular velocity (b) at the end of each simulation for selected models from the si 1 model series, evolved in spherical 
symmetry with SESAME. The models are initially forced into rotation via the rotation law given in eq. with a choice of Qq = 27rPo"' and A as specified in the 
model name. Model graphs plotted in solid lines correspond to spherically symmetric evolution without the inclusion of viscous effects (:1D suffix in TablesfTI 
andl2l. Model graphs plotted in dotted lines correspond to spherically-symmetric evolution with viscous dissipation and angular momentum redistribution (:lDv 
suffix in Tables fTl andEl . For models si 1A1000P1.25 (black) and sllAlOOOPS.OO (blue), no :lDv evolution is carried out. Note the approximate solid-body 
rotation of all : ID models out to '^10 km. Also note the factor of 3 to 4 larger region of solid-body rotation realized by viscous dissipation in all :lDv models. 




Fig. 6. — (a) Final moment-of-inertia-weighted mean peiiod Pm for mateiial interior to a given mass coordinate M and (b) mass M(R) contained inside a 
radius R. Both quantities are plotted for a subset of the sll:lD and :lDv model series at the end of each calculation. The early and strong explosions of models 
sllA1000P2.00:lDv and sllA1000P3.00:lDv lead to significantly smaller PNS masses. This is apparent in plot (b), but also in plot (a), where the lower "mass 
cut" in Pm is obvious. Models sllAl()()0P1.25:lD and s11AcxjP10.47:1D, and their :lDv variants, also explode, but more weakly and at later times. This is 
why their profiles do not deviate as much from the models that do not explode. Model sllA1000P1.25:lD is the fastest model that we consider. Due to strong 
centrifugal forces, si lAl()00P1.25:lD's protoneutron star is rather extended (plot b). Also note that the inclusion of dissipative angular momentum redistribution 
leads to rigid rotation throughout the PNS for all :lDv models. 
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Fig. 7. — Logarithmic color map of tlie density for a variety of 2D models at 200 ms after core bounce. Sliown is the inner 400 km on a side. This figure 
nicely demonstrates the pronounced rotational flattening of the protoneutron star in fast models (si 1A1000P2.34 and EISA). With increasing precollapse iron 
core period, the PNS becomes less and less rotationally flattened. In models sllAlOOOPlO.47 and ml5b6 the flattening has almost disappeared. The lighter 
cores (sll models) show considerably stronger convective structure — with the exception of the purely hydrodynamic variant of si 1A1000P10.47 — than the 
more massive cores (E15, ml5b6) at 200 ms after bounce. The slight "dents" and "spikes" along the rotation axis are artefacts of the imperfect symmetry axis 
treatment. Note that the pure hydro variant of sllAlOOOPlO.47 was evolved with an updated version of the VULCAN/2D with improved axis treatment. 
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Fig . 8. — Mass contained inside a spliere of radius R for selected two-dimensional models ~200 ms after core bounce. The "s" models IWooslev & WeaveJ 
119951) with Po = 2.34 s and the fast E15A and E20A models (see model parameter details in Table ITT deviate significantly from spherical symmetry. Their 
spherical mass profiles should, hence, be taken with caution. Solid lines represent mass profiles of MGFLD models. Dashed lines coiTespond to models which 
are evolved purely hydrodynamically and experience prompt explosions that blast away the outer core. 
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Fig. 9. — Equatorial profiles of the angular velocity (a) and the period (b) at 200 ms after bounce for each model of our sll:2D model series. The two- 
dimensional PNSs exhibit a low degree of differential rotation inside ~20 km and then become strongly differentially rotating further out. All sllA1000:2D 
models exhibit the same qualitative features. Model sllooPlO.47 (initial solid-body rotation), however, differs from its A = 1000 km counterpart. Its PNS is 
appreciably less compact and the region of quasi-solid-body rotation is more extended. The dashed lines correspond to models that are evolved adiabatically. 
These models promptly explode and do not undergo a phase of strong convection. Hence, their equatorial profiles show less structure. The profiles are plotted 
from 3 km on outwai'ds. Due to the imperfect axis treatment in the version of VULCAN/2D that is used for these calculations, we do not highlight the angular 
velocity and period data inside ~3 km. 
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Fig. 10. — Logarithmic color map of the angular velocity fl for models of our si 1A1000:2D model series evolved in axisymmetry using VULCAN/2D with 
Po = 2.34 s, 4.69 s, 10.47 s, and 41.69 s, corresponding to f!o = 2.68 rad s"', 1.34 rad s"', 0.60 rad s"', 0.15 rad s"'. Shown are the inner 400 km on a side. 
Fluid velocity vectors are superposed. The vector data are downsampled at radii larger than ~200 km. The panels show the angular velocity distribution at about 
200 ms after bounce in each model. The fastest model, si 1A1000P2.34:2D, has a large centrifugally supported region of high angular-velocity material. This 
region grows increasingly smaller with increasing initial iron core period. The color map has been artificially limited to a maximum of lO' rad s"' to cut out 
unphysically large angular velocities in the zones very close to the rotation axis. 
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Fig. 1 1 . — Postbounce evolution of the moment of inertia-weighted mean period Ppi2 ■ Solid-line graphs depict standard models of the si 1A1000:2D model 
series. Plotted with the dashed lines are the Ppi2 of models that were evolved purely hydrodynamically. As expected, all Ppi2 exhibit a local minimum at bounce 
when the inner core is most compact. Ppi2 increases slightly after bounce, but ultimately decreases as more material settles onto the contracting PNS. With 
increasing initial iron core spin period, Ppi2 increases systematically. The purely hydrodynamic models' PNSs cannot contract. Hence, their Pp[2 s stay constant 
or increase slightly as lower angular momentum material settles above the density threshold. The fast pure hydro model (Pq = 2.34 s) exhibits strong postbounce 
oscillations like those described in the study of iOtt et all 120041) . Such oscillations are critically damped in simulations including adequate neutrino treatment. 



27 



1 ^= T/|W| for the sll-2D Model Series 

» :SD (solid) and :2Dh (dashed) 

I I 




100 300 300 

t - (ms) 



Fig. 12. — Evolution of the rotation parameter /3 = T/|W| at postbounce times for a subset of the sllA1000:2D model seiies. (3 increases dramatically at 
bounce and increases monotonically, but linearly after bounce in all MGFLD models, because of angular momentum conservation (and hence spin-up) of material 
that settles down onto the PNS. The slope of /3 increases with decreasing precollapse spin period. Despite having the same initial spin (Pq = 10.47 s), model 
si 1AcxdP10.47:2D, which is initially in solid-body rotation, reaches higher postbounce /3s than model si 1A1000P10.47:2D. This is clearly related to the greater 
total angular momentum in the former Model si lA1000P2.34-2Dh's (3 evolution reflects this model's large-scale coherent postbounce oscillations. Both :2Dh 
models evolve to constant (3. both because they explode and because their PNSs cannot contract and spin up. 
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Fig. 13.— Logarithmic color map of the angular velocity distribution in models si 1A1000P2.34:2D, si lA1000P2.34:2Dh, sl5A1000P2.34:2D, and 
s20A1000P2.34:2D. Sh own are the inner 400 km on a side at 200 ms after bounce (150 ms for si lA1000P2.34:2Dh). Velocity vectors are supeiposed in 
the same way as in Fig. 1101 Model si lA1000P2.34:2Dh is evolved purely hydrodynamically. It explodes and leaves behind a less compact configuration with 
smaller central f2 than the corresponding model with MGFLD. With increasing progenitor iron core mass (TablelTT, the magnitude of the final angular velocity 
increases and its distribution becomes broader. 
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Fig. 14. — Logarit hmic color map of the a ngular veloc i ty of m odels E15A:2D, E15A:2Ds ("shellular rotation"), E20A:2D and ml5b6:2D which use the 
rotating progenitors of IHeger et alj 1200(1) and lHeger et alj 12003) . Shown are the inner 400 km on a side at 200 ms after core bounce. Velocity vectors are 
superposed and downsampled at radii greater than about 200 km. The :2Ds variant of E15A is initially set up to rotate with constant Q on spherical shells. Its 
region of high angular velocity is more compact than in its :2D counteipart, which is initially set up to rotate with constant angular velocity on cylindrical shells. 
Interestingly, the innermost regions of E15A:2D and E15A:2Ds are quite similar Model E20A:2D has a more massive core than model E15A:2D, but slower 
initial rotation. Its central angular momentum distribution is very similar' to the one seen in E15A:2D. Model ml5b6:2D has a high initial period of '^32 s. 




Fig. 15. — Equatorial profiles of the angular velo city at 200 ms aft er bounce for a set of models with similar initial iron core periods (Pq) and different 
progenitor mass. The slow "magnetic" model ml5b6 IHeger et all2005j) stands out and has been included in (a) for comparison, (a) shows f7 on a logarithmic 
scale vs. radius and (b) displays f2 on a linear scale for the innermost 100 km of each model. On a logarithmic scale it is hardly possible to distinguish between 
individual fast models. With the exception of si 1 and sl5 that have very similar progenitor stmcture, more massive progenitors lead to slightly higher PNS central 
SI. The model E15A:2Ds (dashed-magenta graph) is setup in "shellular" rotation. This, however, does not lead to a con siderably different equatorial Q profile at 
postbounce times than that of model E15A:2D. Differences in off-equatorial regions are, however, present (see Fig. 1141 . 
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Fig. 16. — Postbounce evolution of the moment of inertia-weighted mean period Ppi2 for models with differing initial iron core spin and progenitor mass. Ppi2 
strongly depends on the initial iron core spin and weakly on progenitor mass and structure. Interestingly, the choice of rotation law ("cylindrical" vs. "shellular" 
— see dashed plot for E15A:2Ds) has little effect on the mean period. The small discontinuities in some of the graphs at eai'ly postbounce times are due to the 
finite resolution of our numerical grid and occur when an entire mass shell first exceeds the density threshold (p > lO'^ g cm"^). 




Fig. 17. — Evolution of the rotation parameter f3 = T/|W|. Here, we compare a representative set of models from our two-dimensional model series. Due to 
the similarity of their progenitor structures, si 1A1000:2D and sl5A1000:2D have almost identical /3 evolutions for the same initial iron core spin period (red and 
green graphs). /3 increases linearly after bounce in all lWoo slev & Weaver 1 1995) models. Models E15A:2D, E15A:2Ds, and E20A:2D have considerably more 
complicated rotational stnictures (Fig.Q. Hence, an appreciably different /3 evolution can be expected (T cx: C^). Model s20A1000P2.34 shows the steepest 
(3 slope and is likely to reach the classical limit for secular rotational instability (/3 > 14 %). However, a low-T/|W| instability might set in at much lower /3 
lOtt et aH200l . 
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Fig. 18. — Contrasting one-dimensional to two-dimensional results: In panel (a), the equatorial angular velocity profiles at 200 ms after bounce of three 
models, namely sllA1000P2.34 (red), sllAlOOOPlO.47 (blue), and E15A (green), each in spherical symmetry (solid lines) and axisymmetry (dashed lines), 
are compared. The profiles are plotted from 5 km on outwards to excise axis artefacts. For E15A:2D, we also plot the "shellular" variant. Panel (b) shows the 
evolution of the moment-of-inertia weighted mean period Ppi2 for the selected set of models. In Panel (c), we display the evolution of the rotation parameter f3 
and panel (d) depicts the mass-(spherical)radius relationship at 200 ms after bounce. 
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Fig. 19. — Mapping from the initial central iron core spin period to the mean PNS period at ~450 ms after bounce for the one-dimensional and ~200 ms 
after bounce for the two-dimensional models. We choose these times for consistency among models from each model set. Most two-dimensional models are 
evolved beyond ~200 ms after bounce while most one-dimensional models are evolved beyond ~450 ms after bounce (see Tablesl2landl3l. On the abscissa, we 
plot the initial central period in seconds and on the ordinate we plot the moment-of-inertia-weighted mean period of the PNS. Models evolved with radiation- 
hydrodynamics in spherical symmetry are marked by sohd boxes. Spherically-symmetric models that include viscosity are marked by open rhombi. Open stars 
symbolize axisymmetric models. In addition, the model icons are color-coded according to their progenitor type, sll models are represented by red, sl5 by 
green, and s20 by blue. Models from the recent studies of Heger et al. 1 2000) and Heaer et al. ( 2005|) are marked black. For the models considered in this study, 
rotational effects are not dominant and the mapping between initial and final rotation period is approximately hnear. 
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Fig. 20. — Equatorial specific angular momentum (j) profiles at 200 ms after core bounce plotted for a representative set of models evolved in axisymmetry 
with VULCAN/2D. Drawn in solid lines are graphs of models evolved with MGFLD and dashed lines represent models that are evolved purely hydrodynamically 
and undergo prompt explosions. At 200 ms after b ounce, the PNS i s already relat ively compact and most of the mass on the grid is located interior to ~50 km 
(Fig.|8j. Models that make use of the E15A, E20A IHeg er et al.'2000), and ml5b6 IHeger et all20O5l) presupernova models show drops in their j profiles at large 
radii which correspond to the discontinuities in their precollapse rotational profiles (Fig.lll. 



